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VELOCITY  PROFILE  CHARACTERIZATION 
FOR  THE  5 -CM  AGENT  FATE  WIND  TUNNELS 


1.  INTRODUCTION 

The  Agent  Fate  Program  is  sponsored  by  the  Defense  Threat  Reduction  Agency 
as  a  Defense  Technical  Objective:  CB.42  -  Environmental  Fate  of  Agents.1  The  objective  of  the 
program  is  to  measure  and  understand  the  physico-chemical  processes  for  Chemical  Warfare 
agents  (CWAs)  on  surfaces  in  order  to  predict  their  persistence  and  fate  in  operational  scenarios 
via  agent  fate  models.  Of  prime  interest  are  the  amount  of  liquid  volatized  and  the  residual  agent 
remaining  in  the  substrates  under  a  set  of  controlled  environmental  conditions  defined  by  wind 
speed,  temperature,  and  relative  humidity.  The  resulting  experimental  data  are  then  used  to 
develop  models  to  predict  threat  persistence  for  operational  scenarios. 

One  of  the  major  elements  of  this  program  consists  of  experiments  performed  in  a 
series  of  5-cm  Wind  Tunnels.*'  These  facilities,  which  fit  in  a  standard  chemical  fume  hood,  are 
used  to  measure  the  release  and  retention  of  CWAs  from  various  materials  under  simulated 
environmental  conditions.  A  key  feature  of  this  wind  tunnel  is  that  it  creates  specified  vertical 
velocity  profiles  representing  the  lower  portion  of  the  velocity  profile  produced  by  a  wind- 
induced,  atmospheric  boundary  layer.  The  attention  to  duplicating  the  detailed  velocity  profile 
over  the  drop/substrate  was  one  of  the  main  features  that  set  this  experimental  arrangement  apart 
from  previous  agent  fate  type  studies.  It  is  important  that  these  velocity  profiles  be  characterized 
in  order  to  validate  and  interpret  the  subsequent  experimental  data  obtained  in  the  5-cm  Wind 
Tunnel,  as  well  as,  to  allow  comparison  with  the  velocity  profiles  present  in  other  different  sized 
wind  tunnels  supporting  the  Agent  Fate  program.  Characterization,  as  used  here,  refers  to  the 
measurement  of  the  velocity  profiles  over  the  length  and  width  of  the  test  section  under  various 
test  conditions.  This  report  summarizes  the  instrumentation  employed,  the  characterization 
procedure  and  resulting  velocity  profiles  measured  in  the  Version  3  and  Version  3,  Mod  1 
configurations  of  the  5-cm  Wind  Tunnel.  These  tunnel  configurations  were  used  for  the  HD  on 
glass  (tunnel  validation)  and  the  HD  on  sand  test  phases  of  the  Agent  Fate  Program.5 


2.  BACKGROUND 

The  Agent  Fate  Wind  Tunnel  studies  are  intended  to  provide  data  to  the  modelers 
on  the  fate  sequence  of  agent  drops  on  different  substrates  under  the  combined  influence  of  wind, 
temperature,  and  humidity.  Figure  1  describes  the  overall  fate  of  agent  deposited  on  a  substrate 
and  illustrates  the  sequence  of  events  in  terms  of  the  agent  vapor  and  agent  mass  being  measured 
in  the  wind  tunnels.  The  figure  does  not  contain  actual  experimental  data,  but  is  intended  to 
illustrate  the  types  of  data  obtained  in  the  wind  tunnel  test. 
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Figure  1.  Agent  Fate  Sequence 

For  the  case  shown  in  the  figure,  the  substrate  is  both  sorptive  and  reactive.  At 
the  start  of  the  test,  a  drop  or  drops  of  agent  are  deposited  on  a  substrate  specimen.  A  portion  of 
the  agent  forms  a  sessile  drop  on  the  surface  while  the  remainder  is  sorbed  into  the  substrate. 

The  substrate  sample  is  then  inserted  into  the  floor  of  the  wind  tunnel  test  section.  The  wind 
tunnel  air  flow  is  then  established  for  the  desired  specific  velocity,  temperature,  and  relative 
humidity  conditions.  The  chemical  volatilizes  into  the  air  passing  over  the  substrate  with  the 
resulting  vapor  carried  away  by  the  air  stream.  As  time  passes,  the  sessile  drop  on  the  surface 
completely  volatilizes  leaving  only  the  chemical  sorbed  in  the  substrate  which  also  begins  to 
desorb.  This  phase  exhibits  a  markedly  lower  volatilization  rate  than  was  the  case  for  the  sessile 
drop.  Testing  continues  until  no  vapor  appears  to  be  released  from  the  substrate.  This  may 
range  from  minutes  to  weeks  depending  on  the  chemical,  substrate  material,  drop  size,  velocity, 
temperature,  and  relative  humidity. 


Throughout  this  operation,  the  agent  vapor  and  air  mixture  are  analyzed  by  means 
of  an  on-line,  Thermal  Desorption  Gas  Chromatography  -  Mass  Spectrometer  (TD-GC-MS) 
system  and/or  by  sorbent  tubes  to  determine  the  agent  vapor  mass  concentration  released  as  a 
function  of  time.  The  amount  of  agent  released  into  the  air  represents  important  data  for  the 
model  development  because  it  represents  the  initial  conditions  for  transport  and  diffusion 
models.  These  models  predict  the  movement  of  the  vapor  cloud  under  the  influence  of  wind 
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over  the  terrain.  The  TD-GC-MS  can  also  be  used  to  determine  whether  any  of  the  agent 
decomposes  in  the  air.  After  the  test,  the  substrate  can  be  analyzed  for  any  residual  chemical  or 
chemical  decomposition  products  by  means  of  various  analytical  methods  including 
spectroscopy  (e.g.,  in  situ  NMR),  extraction  followed  by  chromatographic  analysis  or  other 
methods.  In  this  way,  the  fate  of  the  chemical  from  each  of  these  mechanisms  can  be  determined 
for  the  experimental  conditions  evaluated.  Any  chemical  remaining  in  the  substrate  is  also  a 
major  concern  to  model  development  because  it  may  represent  a  contact  hazard  or  be  released  at 
a  future  time  due  to  surface  displacement,  rain,  or  other  effects. 


3.  5-CM  WIND  TUNNEL  CONFIGURATIONS 

3.1  Tunnel  Version 


The  Version  3  Wind  Tunnel  configuration,  shown  in  Figure  2,  is  an  open  circuit 
arrangement  with  a  “Twisted  S”  shape  layout.  Figures  3  and  4  contain  schematics  of  the  wind 
tunnel  denoting  its  major  components  consisting  of  a  transition  cone,  turning  vane  section,  fetch, 
test  section,  mixing  section,  static  mixer,  sampling  section  and  support  frame.  The  tunnel  has  an 
overall  width  of  93.2  cm,  an  overall  height  of  56.7  cm,  and  an  overall  depth  of  55  cm.  Figure  5 
depicts  a  cut-away  view  revealing  the  internal  elements  in  the  turning  vane  section,  fetch,  test 
section  and  exhaust  section. 

Except  for  the  nickel  coated,  plastic  transition  cone,  all  tunnel  components  are 
fabricated  from  stainless  steel.  The  fetch,  mixing  section  and  exhaust  section  are  made  from 
5-cm  external  dimension,  square  cross-section,  stainless  steel  tubing.  While  the  inside 
dimension  of  the  tubing  is  4.85  cm,  the  tunnel  is  referred  to  as  a  5-cm  Wind  Tunnel.  The  static 
mixer  and  test  section  are  constructed  from  stainless  steel  sheet  and  plate,  respectively.  Many  of 
the  stainless  steel  components  included  welded  flanges,  allowing  them  to  be  connected  by  screws 
with  silicon  rubber,  chemical  resistant,  Kalrez  Rr  gasket  material  in  between. 
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Figure  2.  Version  3  (Production)  Wind  Tunnel 
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Figure  3.  Version  3  (Production)  Tunnel  Schematic:  Top  View 
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Figure  4.  Version  3  (Production)  Tunnel  Schematic:  Front  View 
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Deleted  for  Clarit> 


Figure  5.  Cut-Away  View  Showing  Internal  Elements 


A  Miller-Nelson  Environmental  Control  Unit  (ECU)  provided  the  tunnel  air  flow 
at  a  prescribed  temperature  and  relative  humidity.  The  unit  was  located  in  front  of  the  fume 
hood  and  supplied  air  through  a  1  cm  internal  diameter  tube  into  the  transition  cone.  The  turning 
vane  section,  shown  in  Figure  6,  contains  three  constant  radius  vanes  that  direct  the  air  flow 
through  a  90°  bend  into  the  entrance  to  the  fetch.  This  arrangement  resulted  in  the  “twisted  S" 
layout  of  the  tunnel. 

The  floor  of  the  fetch  entrance  contains  turbulence  generators  that  alter  the  air 
flow  to  create  the  desired  velocity  profiles  in  the  test  section.  These  are  composed  of  an  ordered 
array  of  23  Teflon  cubes  (4.8  mm  on  a  side)  glued  to  a  flat  stainless  steel  strip  glued  to  the  floor 
of  the  fetch  as  shown  in  Figure  7.  Considerable  development  effort  went  into  evolving  this 
turbulence  generator  array  as  discussed  in  Section  4.2. 
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Figure  6.  Turning  Vane  Section 


Figure  7.  Cubical  Turbulence  Generators:  Top  View 
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The  test  section,  shown  in  Figure  8,  is  fabricated  from  machined  stainless  steel 
plates  attached  by  screws.  Linear  ethylene  propylene  rubber  gaskets  (fabricated  from  cut-up 
O-rings)  provided  a  seal  between  the  plates  interfaces.  The  test  section  ceiling  contains  a  total  of 
seven  instrument  ports.  These  include  an  array  of  five  ports  for  the  hot  wire  anemometer  probe: 
a  central  location,  an  upstream  and  a  downstream  location,  and  a  left  and  a  right  lateral  location. 
All  five  positions  are  used  for  the  velocity  characterization  surveys.  The  port  location 
downstream  of  the  center  position  was  used  for  the  reference  velocity  probe  during  the  agent  fate 
experiments.  The  test  section  temperature  is  monitored  through  the  port  located  most 
downstream  from  the  center  port  position.  A  4.85  cm2  flat  side  window  is  located  on  the  side  of 
the  test  section.  It  is  flush  with  the  test  section  floor  to  provide  an  unobstructed  side  view  of  the 
sessile  drop  of  agent  on  the  substrate  surface. 


Mixing  Section  Instrumentation  Fetch  Section 

Ports 


Air  Flow 


Window  Piston 

Housing 


Substrate  Test  Section 

Surface 


Figure  8.  Test  Section 


Prior  to  a  test,  the  substrate  is  placed  on  the  upper  surface  of  a  cylindrical  piston 
located  in  the  floor  of  the  test  section.  The  piston,  which  is  composed  of  a  Teflon  center  section 
and  a  stainless  steel  outer  body,  is  located  in  the  center  of  the  test  section  floor.  The  outer  body 
is  threaded  and  is  inserted  into  a  threaded  hole  in  the  test  section  floor  until  the  top  surface  of  the 
substrate  is  flush  with  the  test  section  floor.  The  top  surface  of  the  piston  measures  3.8  cm  in 
diameter  and  can  contain  either  a  circular  or  square  cavity  to  hold  the  substrate  sample.  The 
lower  portion  of  the  piston’s  outer  body  includes  a  knurled  finish  to  facilitate  installation  and 
removal  of  the  piston.  A  rubber  “O”  ring  is  located  in  a  circumferential  slot  on  the  outer  surface 
of  the  piston  as  a  seal  between  the  piston  and  the  test  section  floor. 
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The  static  mixer  is  a  hollow,  cubical  box  having  internal  dimensions  of  10  cm  on 
a  side,  considerably  larger  than  either  the  inlet  or  exit  ducts.  The  air  flow  from  the  mixing 
section  enters  the  mixing  box  where  it  circulates  and  mixes  before  exiting  at  a  90°  angle  into  the 
exhaust  duct.  The  vapor  sampling  probe  consists  of  a  stainless  steel  tube  welded  to  the  exhaust 
and  enters  from  the  side,  curving  through  a  90°  bend  to  face  directly  into  the  air  flow.  The 
flexible  Teflon  inlet  tubing  leading  to  the  sampling  instrumentation  is  inserted  into  the  sampling 
tube  until  the  leading  edges  of  the  flexible  tube  is  flush  with  the  end  of  the  stainless  steel  tube. 
This  placed  the  sampling  inlet  5  cm  upstream  from  the  exit  of  the  exhaust  section,  on  the 
centerline  of  the  exhaust  duct  and  oriented  parallel  to  the  air  flow. 

The  primary  vapor  sampling  instrument  is  the  Versatile  Tube  Sampler  (VTS).' 
The  VTS  consists  of  a  series  of  30  vapor  sorbing  Depot  Area  Air  Monitoring  System  (DA  AMS) 
tubes  that  sequentially  collect  the  vapor/air  mixture  at  the  tunnel  sample  port.  After  the  test,  the 
content  of  the  tubes  are  analyzed  providing  a  time  history  of  vapor  concentration.  In  order  to 
maintain  constant  temperature  conditions  throughout  the  wind  tunnel,  its  entire  structure  is 
enclosed  in  an  insulation/temperature  controlled  blanket.' 
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Figure  9.  Version  3  Tunnel  and  Associated  Instrumentation5 
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Details  of  the  Version  3  tunnel  and  its  associated  instrumentation  are  depicted  in 
Figure  9.  Figure  10  illustrates  the  tunnel  installed  in  the  fume  hood  along  with  a  person  as  a  size 
reference. 


5-cni  Wind  Funnel 


Fume  Hood 
Window 


Figure  10.  Version  3  Tunnel  Installed  in  Fume  Hood 

Depending  on  the  desired  test  temperature  conditions,  the  wind  tunnel  side  walls 
are  either  heated  or  cooled.  For  the  medium  and  high  test  temperature,  both  which  are  above 
room  temperature,  the  tunnel  is  wrapped  with  electric  heating  tape.  The  configuration  of  a 
tunnel  for  low  temperature  tests  require  wrapping  the  upstream  sections  up  to  and  including  the 
test  section  with  copper  tubing  through  which  cold  water  was  circulated.  Thermocouples 
monitored  and  controlled  the  tunnel  wall  temperature  which  had  to  remain  constant  throughout 
the  entire  test.  Because  of  the  special  arrangements  required  to  achieve  these  two  ranges  of 
temperatures,  certain  tunnels  were  dedicated  for  medium  and  high  temperature  tests  and  others 
for  cold  temperature  tests.' 

3.2  Tunnel  Version  3  Mod  1 


The  Version  3  Mod  1  tunnels  are  intended  for  installation  in  the  large  capacity 
chemical  fume  hoods.  A  new  support  structure  was  designed  allowing  a  group  of  four  Version  3 
tunnels  to  be  installed  in  a  single  hood.  This  arrangement  required  the  humidity  and  flow  control 
units  to  be  located  to  the  side  of  the  hood  permitting  the  hood  sash  to  be  closed  during  testing 
operations.  In  order  to  accomplish  this,  the  transition  cone  was  placed  in  a  vertical  orientation 
giving  the  tunnel  a  “U”  shaped  layout  as  shown  in  Figure  1 1 . 
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Figure  1 1.  Version  3  Mod  1  Wind  Tunnel 


The  horizontal  orientation  of  the  turning  vane  section  of  the  Version  3  tunnel 
configuration  resulted  in  any  flow  disturbances  due  to  the  vanes  to  act  in  a  horizontal  plane. 
While  this  influenced  the  variance  of  the  profiles  in  a  lateral  sense,  it  had  little  effect  on  the  more 
critical  shape  of  the  vertical  velocity  profiles.  However,  the  vertical  orientation  of  the  turning 
vane  section  of  the  Version  3  Mod  1  tunnel  configuration  resulted  in  any  flow  disturbances  due 
to  the  vanes  to  act  in  a  vertical  plane  and  consequently  had  a  greater  influence  on  the  shape  of 
the  vertical  velocity  profiles.  Based  on  velocity  profile  measurements,  it  was  determined  that  the 
Version  3,  Mod  1  tunnel  required  an  inlet  screen  to  eliminate  air  flow  perturbations  produced  by 
the  turning  vanes.  A  24  x  24  mesh  per  inch,  0.014  in.  diameter  wire  screen  was  located  at  the 
exit  of  the  turning  vane  section  (the  entrance  of  the  fetch)  as  shown  in  Figure  12. 
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Figure  12.  Fine  Mesh  for  Vertical  Oriented  Turning  Vane  Section 

Another  modification  incorporated  into  the  Version  3  Mod  1  tunnel  involved 
replacing  the  screw-in  piston  with  a  snap-in  piston.  This  improved  design  allowed  a  more  rapid 
insertion  and  removal  from  the  tunnel  than  the  previous  screw-in  piston  design.  This  greatly 
reduced  the  time  for  the  temperature  to  recover  and  for  the  start  of  the  measurement  of 
evaporation  data  during  the  crucial  early  portion  of  the  test. 

In  preparation  for  tests  with  VX,  a  low  volatility  agent,  all  Version  3  Mod  I 
tunnels  were  electro-polished  and  SulfInertK  coated.  The  SulfInertK  treatment  provided  an  order 
of  magnitude  increase  in  resistance  to  absorption  of  VX  vapor  compared  with  the  Silco-Steef 
process.  The  Version  3  Mod  1  tunnel  is  being  used  for  the  major  portion  of  the  Agent  Fate 
production  testing. 

A  total  of  six  Version  3  Wind  Tunnels  were  originally  fabricated  (designated  3  A 
through  3F)  and  three  of  these  were  used  for  the  5-cm  Wind  Tunnel  validation  tests  with  HD  on 
glass  and  HD  on  sand.  A  total  of  six  Version  3  Mod  1  tunnels  were  fabricated  and  designated 
3G  through  3L.  Table  1  summarizes  the  current  disposition  of  the  twelve  5-cm  Wind  Tunnels 
and  the  status  of  their  velocity  characterization.  The  term  “Complete"  in  the  Velocity 
Characterization  column  means  that  vertical  velocity  profiles  were  measured  for  all  three  wind 
speeds  at  all  five  horizontal  locations  relative  to  the  center  of  the  substrate  piston  (i.e.,  center, 
upstream,  downstream,  left  and  right).  The  term  “Abbreviated”  means  that  a  subset  of  velocities 
and  positions  were  characterized.  In  the  future.  Version  3  Wind  Tunnels  may  be  upgraded  to 
Version  3  Mod  1  Tunnels. 
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Table  1.  Disposition  and  Use  of  5-cm  Wind  Tunnels 
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4. 


VELOCITY  PROFILES 


4. 1  Specified  Operational  Velocity  Profiles 

The  approach  adopted  by  the  Agent  Fate  program  was  to  define  a  set  of  nominal 
vertical  velocity  profiles,  based  on  boundary  layer  theory  and  the  “operational'’  wind  conditions 
of  interest  to  the  Agent  Fate  model  developers.6  While  somewhat  arbitrary  for  simplified 
conditions,  these  profiles  represent  a  practical  bridge  between  the  variable  turbulent,  unsteady, 
and  omni-directional  air  flow  in  nature  and  the  turbulent,  steady,  and  parallel  flow  of  the  wind 
tunnel. 


These  velocity  profiles  were  developed  using  boundary  layer  theory  and  are  based 
on  three  values  of  Friction  Velocities.  Friction  Velocity  relates  the  shear  stress  produced  by  a 
boundary  layer  on  the  surface  of  the  ground  and  is  defined  in  Figure  13. 


I T  -  Friction  Velocity  =  \Txx/p 

txx  =  Surface  shear  stress  -  p  du/d>  |%  0 

p  -  Viscosity  of  air 

du/dyl,  H  -  velocity  gradient  at  \  0 

p  -  Air  density 

u  =  Velocity 

y  -  Height 


>=0 


Figure  13.  Velocity  Profile  Terminology 

A  set  of  three  “operational”  boundary  layer  velocity  profiles  were  defined  by  the 
model  developers  as  standard  conditions  for  the  wind  tunnel  studies.  Figure  14  contains  the 
original  graph  provided  by  the  model  developers  and  defines  the  near  surface  velocity  profiles  of 
an  atmospheric  boundary  layer  on  a  flat  surface  with  some  degree  of  roughness  for  wind  speeds 
of  0.5,  3.0,  and  6.0  m/s  at  a  2  m  height  above  the  ground.  Wind  velocity  conditions  at  2  m  were 
selected  based  on  operational  conditions  intended  to  be  included  in  the  final  Agent  Fate  models. 
The  Friction  Velocities  selected  for  the  Agent  Fate  experiments  correspond  to  relatively  smooth 
surfaces  (glass,  metal,  sand,  soil,  concrete  and  asphalt).  The  choice  of  these  friction  velocities 
defined  their  respective  velocity  profiles. 
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Figure  14.  Operational  Velocity  Profiles 


Figure  15  shows  the  three  operational  velocity  profiles  in  non-logarithm  units  and 
provides  a  better  physical  “feel”  for  the  shape  of  the  velocity  profile.  While  the  operational 
velocity  profile  is  defined  for  a  wind  velocity  at  a  height  of  2  m,  the  5-cm  Wind  Tunnel  only 
needs  to  duplicate  the  velocity  profile  up  to  a  height  of  2  cm  as  shown  in  Figure  16.  The 
velocities  at  a  height  of  2  cm  in  the  5-cm  Wind  Tunnel  of  0.22,  1.7  and  3.6  m/s  correspond  to  the 
0.5,  3.0  and  6.0  m/s  test  matrix  velocities,  respectively.  Hot  wire  anemometer  measurements  in 
the  5-cm  Wind  Tunnel  test  section  had  demonstrated  that  the  stipulated  operational  vertical 
velocity  profiles  were  created  at  these  reference  velocities.  Thus,  setting  the  2-cm  reference 
velocity  at  a  specific  test  matrix  value  was  sufficient  to  establish  the  associated  operational 
vertical  velocity  profile.  The  reference  height  was  later  lowered  to  1  cm  for  the  production  tests. 
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Figure  15.  Operational  Velocity  Profiles  in  Physical  Units 
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Figure  16.  Lower  Portion  of  Operational  Boundary  Layer  Profiles 


26 


The  enlarged  plot  of  the  profiles  contained  in  Figure  16  shows  the  characteristic 
(and  unrealistic)  “dog  leg”  shape  of  the  low  operational  velocity  profile  in  the  vicinity  of  the 
floor.  This  is  due  to  the  manner  in  which  the  transition  between  the  laminar  sub-layer  and  the 
exponential  portions  of  the  theoretical  profile  was  established  when  the  profiles  were  defined. 
The  measured  (and  realistic)  shape  exhibits  a  smoother  curve.  A  possible  explanation  for  this 
unrealistic  shape  may  be  due  to  the  application  of  turbulent  boundary  layer  equations  to  the  low 
velocity  case,  where  turbulent  flow  is  not  possible. 

4.2  Formulation  of  5-cm  Wind  Tunnel  Velocity  Profiles 


The  5-cm  Wind  Tunnel  utilizes  a  series  of  turbulence  generators  located  upstream 
of  the  test  section  to  produce  the  specified  operational  vertical  velocity  profiles  over  the 
agent/substrate.  These  turbulence  generators  consisted  of  an  array  of  cubical  blocks  placed  along 
the  floor  at  the  beginning  of  the  fetch.  The  turbulence  created  by  these  blocks  produce  vertical 
velocity  profiles  that  closely  match  the  three  operational  profiles.  The  location  of  the  turbulence 
generators  relative  to  the  test  section  is  described  in  Figure  17  and  details  of  the  turbulence 
generator  array  are  shown  in  Figure  1 8. 


5.  5-CM  WIND  TUNNEL  VELOCITY  CHARACTERIZATION 

5. 1  Overview  of  Velocity  Characterization 


“Characterization"  refers  to  the  detailed  measurement  of  the  vertical  velocity 
profiles  in  the  test  section  of  the  5-cm  Wind  Tunnel  and  relating  these  profiles  to  a  reference 
velocity  that  could  be  easily  measured  in  the  tunnel  during  routine  production  tests.  The 
characterization  process  was  complicated  by  the  small  size  of  the  test  section,  the  relatively  low 
velocities  and  the  large  range  of  temperature  and  relative  humidity  conditions,  and  the  limited 
number  of  access  ports.  Selection  of  the  velocity  measurement  instrumentation  and  its 
subsequent  use  was  a  critical  element. 
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Figure  17.  Location  of  Turbulence  Generators  Relative  to  Test  Section 
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Figure  18.  Details  of  Turbulence  Generator  Array 
(dimensions  in  cm) 


5.2  Characterization  Instrumentation 

5.2.1  Hot  Wire  Anemometers 

5. 2. 1.1  Overview  Hot  Wire  Anemometry 

Several  techniques  were  considered  to  measure  the  velocity  profiles  including 
pitot-static  tubes,  Particle  Imaging  Velocimeter  (PIV),  Laser  Doppler  Anemometry  (LDA)  and 
hot  wire  anemometry.  Each  had  advantages  and  disadvantages,  but  hot  wire  anemometry 
provided  the  best  approach  for  meeting  the  velocity  profile  measurements  required  for  the  5-cm 
Wind  Tunnel. 


A  Constant  Temperature  Anemometer  (CTA)  operates  by  having  a  current  passed 
through  a  small  wire  and  heating  it  to  200-300  °C.  When  air  flows  over  the  wire,  heat  is  lost  by 
convection  and  the  wire  cools.  The  anemometer  contains  an  electronic  feedback  circuit  that 
increases  the  current/voltage  to  maintain  a  constant  hot  wire  temperature.  The  velocity  of  the  air 
flow  over  the  wire  is  proportional  to  the  voltage  change  required  to  maintain  the  wire 
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temperature.  Calibration  curves  are  generated  relating  air  speed  to  circuit  voltage  for  constant  air 
temperatures.  A  major  difficulty  in  hot  wire  anemometry  is  faced  when  measuring  wind  speed  in 
air  that  is  at  a  different  temperature  from  existing  calibration  curves.  Because  the  CTA  works  on 
heat  transfer,  errors  in  velocity  readings  can  occur  for  air  temperatures  not  accounted  for  by 
either  previous  calibration,  electronic  adjustment  of  the  hot  wire  overheat  ratio,  or  analytical 
corrections  to  the  hot  wire  output  voltage.  Corrections  for  air  temperature  differences  depend  on 
the  type  of  hot  wire  anemometer  employed.  Industrial  type  hot  wire  anemometers  intended  for 
Heating,  Ventilation,  and  Air  Conditioning  (HVAC)  work  typically  well  for  air  temperature 
electronically  using  a  temperature  sensor  placed  near  the  hot  wire  sensor.  Research  grade 
systems  usually  rely  on  a  combination  of  multiple  calibration  curves  at  different  air  temperatures 
and/or  analytical  techniques. 

Both  types  of  hot  wire  anemometers  were  used  with  the  5-cm  Wind  Tunnel.  The 
detailed  characterization  was  performed  with  a  research  grade  anemometer  and  an  industry 
HVAC  grade  anemometer  was  used  to  measure  the  reference  velocity  during  the  actual  testing  in 
the  tunnel. 

5. 2. 1.2  Characterization  Hot  Wire  Anemometer 

Characterization  measurements  were  performed  using  a  TS1  Miniature  Straight 
Hot  Wire  Probe  (Model  1260A).  This  probe  was  calibrated  for  wind  speeds  as  low  as  0.05  m/s. 
A  high  speed  analog  to  digital  converter  was  utilized  along  with  a  dedicated  software  program. 

In  addition  to  velocity,  the  device  also  measures  turbulence  intensity,  Reynolds  stress  and  other 
statistical  quantities.  Multi-component  hot  wires  are  available  which  allow  velocities  to  be 
measured  in  two  and  three  components,  these  are  larger,  more  complicated  to  calibrate  and  use 
and  do  not  offer  as  fine  a  spatial  resolution  as  a  single  component  probe.  Accordingly,  all  of  the 
5-cm  characterization  measurements  were  done  with  a  single-component  hot  wire.  This  was  also 
felt  to  be  sufficient  in  that  the  air  flow  in  the  5-cm  Wind  Tunnel  test  section  is  essentially  uni¬ 
directional.  Figure  19  contains  a  photograph  of  the  basic  TSI  Model  1260.  A  hot  film  probe 
with  its  dimensions  is  shown  in  Figure  20. 
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Figure  19.  TS1  Miniature  Straight  Hot  Wire  Probe  (Model  1260A) 
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Figure  20.  Dimensions  of  TS1  Hot  Wire  Probe  (Model  1260A) 
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The  Model  1260A  hot  wire  anemometer  operates  in  conjunction  with  the  TSI  1FA 
300  controller  shown  in  Figure  21 .  When  air  flows  over  the  wire,  heat  is  lost  by  convection  and 
the  wire  cools.  The  anemometer  contains  an  electronic  feedback  circuit  that  increases  the 
current/voltage  to  maintain  a  constant  hot  wire  temperature.  Through  the  electronics  of  the 
1FA300,  the  temperature  of  the  sensor  is  held  constant  by  increasing  the  current  through  the 
sensor  which  is  outputted  from  the  IFA300  as  a  voltage  difference.  This  voltage  is  supplied  to  an 
analog  to  digital  converter  inside  a  computer.  The  velocity  of  the  air  flow  over  the  wire  is 
proportional  to  the  voltage  change  required  to  maintain  the  wire  temperature.  Specialized 
software  (ThermalPro®)  developed  by  TSI  provides  complete  control  of  the  IFA300  from 
calibration  through  the  actual  measurement  process,  to  velocity  data  analysis.  An  important 
feature  of  the  IFA300  is  its  ability  to  measure  high  frequency  velocity  fluctuations,  which  is 
contributes  greatly  to  the  study  of  turbulent  flow  fields.  The  system  also  incorporates  a 
thermocouple  for  measuring  air  temperature  for  applying  off-calibration  temperature  corrections. 


Figure  21 .  TSI  IFA  300  Thermal  Anemometry  Systems 


It  is  generally  recommended  that  the  probe  be  installed  so  that  the  flow  is  parallel 
to  the  needle  like  hot  wire  supports.  However,  due  to  the  reasons  just  cited,  straight  probes  were 
used  in  characterizations  of  the  ECBC  tunnel  with  the  wire  supports  oriented  perpendicular  to  the 
flow.  This  is  an  acceptable  procedure  for  probes  calibrated  perpendicular  to  the  wind  direction 
as  is  the  case  with  the  5-cm  Wind  Tunnel. 
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5.2. 1.3 


Hot  Wire  Anemometer  Calibration 


A  research  grade  hot  wire  anemometer  requires  accurate  calibration  for  the  range 
of  air  flow  temperatures  for  which  it  is  being  used.  This  calibration  was  performed  using  a  TSI 
Model  8390  bench  top  wind  tunnel  illustrated  in  Figure  22.  The  calibration  wind  tunnel 
incorporates  a  10-cnr  cross-section  test  section  which  can  be  configured  by  adding  or  removing 
supplied  orifice  plates  to  cover  a  manufacturer  specified  wind  speed  range  of 
0.15  to  45.7  m/s.  The  tunnel  air  flow  is  provided  by  a  fan  located  downstream  of  the  test  section 
powered  by  a  variable  speed,  electric  motor.  The  test  section  velocity  had  been  previously 
measured  by  TSI  using  a  laser  doppler  velocimeter  and  these  measurements  are  correlated  to  the 
static  pressure  difference  between  the  settling  chamber  and  test  section.  To  accommodate  the 
expected  low  wind  speeds  near  the  floor  of  the  5-cm  Wind  Tunnel,  the  lower  calibration  wind 
speed  value  was  lowered  to  0.05m/s,  which  was  near  the  lower  limit  of  the  lab's  pressure 
measuring  capabilities.  The  anemometer  being  calibrated  was  inserted  through  the  side  of  the 
calibration  wind  tunnel’s  test  section  so  that  the  tip  of  the  anemometer  was  located  at  the  center 
of  the  test  section. 


Inlet 


Test  Section 


Exhaust 


Figure  22.  TSI  Model  8390  Hot  Wire  Anemometer  Calibration  Wind  Tunnel 


The  detailed  procedure  followed  in  calibrating  the  TSI  Model  1260A  anemometer 
is  described  in  Appendix  I. 


33 


5.2. 1.4 


Reference  Velocity  Hot  Wire  Anemometer 


The  velocity  characterization  was  performed  using  a  research  grade  anemometer. 
This  demonstrated  that  the  wind  tunnel  being  tested  produced  the  required  test  matrix  velocity 
profiles.  During  the  actual  wind  tunnel  tests  of  agents/substrates,  the  desired  velocity  profiles 
were  established  by  setting  the  tunnel  air  flow  through  the  measurement  of  a  reference  velocity 
at  a  specific  height  from  the  test  section  floor.  This  reference  velocity  was  correlated  to  the 
particular  test  matrix  velocity  profile  being  considered.  A  TS1  Air  Velocity  Transducer  Model 
8465  hot  wire  anemometer8  was  employed  to  measure  the  reference  velocity  during  testing 
operations.  A  picture  of  the  Model  8455  is  shown  in  Figure  23. 


Figure  23.  TSI  Air  Velocity  Transducer  Model  8465 

This  device,  while  not  as  sensitive  and  accurate  as  a  research  level  anemometer, 
was  more  robust  and  less  prone  to  damage  or  changes  in  calibration  during  the  long,  sustained 
production  testing  operations  with  the  5-cm  Wind  Tunnel.  It  was  also  thinner  and  was  less 
obtrusive  to  the  air  flow  in  the  passing  over  the  drop/substrate  than  the  research  grade  device. 
The  Model  8465  hot  wire  sensor  is  a  thin  needle  like  wire  that  protrudes  from  the  end  probe  tip. 
The  hot  wire  sensor  is  not  protected  thus  reducing  flow  blockage  which  might  otherwise 
adversely  affect  the  wind  speed  measurement  in  confined  spaces.  A  close-up  of  tip  is  shown  in 
Figure  24. 
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Figure  24.  Schematic  of  Model  8465  Hot  Wire  Probe 


The  Model  8465  is  a  relatively  low  cost,  active  temperature  compensated  device. 
The  small  probe  size  allows  very  accurate  placement  of  the  probe  tip  with  respect  to  the  tunnel 
floor,  while  also  incorporating  some  degree  of  ruggedness  that  research  type  hot  wire  sensors  do 
not.  The  sensor  overheat  ratio  for  the  Model  8465  is  typically  lower  than  those  used  by  research 
hot  wire  systems  thus  reducing  the  chance  of  convection  affecting  the  velocity  readings  near  a 
surface.  The  operational  fluid  temperature  range  of  this  device  is  0  to  60  °C,  with  an  accuracy  of 
±2%  of  the  reading.  However,  this  accuracy  is  only  valid  over  a  temperature  range  of  1 8  to 
28  °C.  Outside  of  this  range  (0  to  18  °C  and  28  to  60  °C)  an  additional  0.2%  per  °C  must  be 
added  to  the  ±2%.  Thus  if  the  air  temperature  was  50  °C,  the  accuracy  would  be  ±6.4%. 

This  anemometer  is  intended  to  make  average  flow  measurements  in  the  center  of 
relatively  large  ducts  and  is  not  suited  for  boundary-layer  measurements.  This  probe  works  best 
when  the  entire  tip  is  immersed  in  the  flow  field  where  the  velocity  is  nearly  constant.  In  a  non- 
uniform  flow  field  such  as  within  a  boundary-layer,  the  probe  would  measure  an  average 
velocity  over  the  length  of  the  10  mm  probe  length.  However,  it  is  adequate  for  reference 
velocity  monitoring  in  the  vicinity  of  the  centerline  of  the  5-cm  Wind  Tunnel  where  the  vertical 
velocity  profile  is  relatively  constant. 

The  tunnel  air  flow  rate  was  established  by  measuring  the  velocity  in  the  test 
section  at  a  reference  height  above  the  test  section  floor  using  a  hot  wire  anemometer  and 
adjusting  the  wind  tunnel  flow  controller  to  achieve  the  test  matrix  reference  velocity  desired. 
Based  on  the  velocity  characterization,  this  provided  the  stipulated  velocity  profile  in  the  test 
section  as  shown  in  Figure  25.  The  velocity  values  at  a  2-cm  height  are  the  source  of  the  0.22, 

1 .7,  and  3.6  m/s  for  the  low,  medium  and  high  operational  velocities,  respectively. 
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Figure  25.  5-cm  Wind  Tunnel  Reference  Velocities 


5.3  Characterization  Arrangement  and  Process 

The  characterization  process  required  a  positioning  device  to  accurately  place  the 
hot  wire  probe  at  specific  heights  above  floor  of  the  test  section.  A  24  in.  commercially 
available  height  gage  (24  in.  Mitutoyo  Series  192)  was  modified  to  position  the  probe  holder 
through  the  instrumentation  ports  in  the  test  section  ceiling.  This  permitted  precise  vertical 
positioning  of  the  hot  wire  sensor  relative  to  the  tunnel  floor.  Figure  26  shows  the  hot  wire 
anemometer  positioning  system  with  the  5-cm  Wind  Tunnel.  A  rubber  balloon  was  used  as  a 
flexible  seal  to  prevent  air  leakage  through  the  port  hole  in  the  test  section  ceiling  and  the  hot 
wire  support  strut.  A  close-up  photo  of  the  anemometer  probe  in  the  tests  section  is  contained  in 
Figure  27. 
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Figure  26.  Hot  Wire  Vertical  Positioning  Apparatus 
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Figure  27.  Hot  Wire  Probe  in  Test  Section 
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A  complete  characterization  consisted  of  measurements  of  the  vertical  profiles  at 
all  five  horizontal  locations  relative  to  the  center  position  defined  as  the  center  of  the  substrate 
piston  (i.e.,  center,  upstream,  downstream,  left  and  right).  The  details  of  the  test  section  and 
measurement  positions  are  shown  in  Figures  28,  29,  and  30. 


Hot  Wire 
Anemometer 
Probe 


10  16 

—  5  08 


Figure  28.  Side  View  of  5-cm  Wind  Tunnel  Test  Section  Showing  Hot  Wire  Anemometer 
Probe  Height  (y)  for  Velocity  Profile  Characterization 
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Figure  29.  Top  View  of  5-cm  Wind  Tunnel  Test  Section  Showing  Hot  Wire  Anemometer 
Probe  Locations  Used  for  Velocity  Profile  Characterization 
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Figure  30.  Positions  for  Vertical  Velocity  Profile  Measurements 


A  typical  characterization  consisted  of  21  velocity  and  turbulence  intensity 
measurements  starting  from  a  location  0.05  in.  (1.25  mm)  above  the  test  section  floor  and 
proceeding  in  increments  of  0.05  in.  ( 1 .25  cm)  to  a  height  of  0.295  in.  (7.5  mm).  From  this 
height  onwards,  the  increment  was  changed  to  0.10  in.  (2.5  mm)  and  the  process  continued  to  a 
height  of  1 .77  in.  (45  mm).  A  detailed  description  of  the  characterization  procedure  is  contained 
in  Appendix  11. 


Figure  31  contains  representative  velocity  measurements  in  the  5-cm  Wind 
Tunnel  at  the  center  position  (A)  compared  to  the  corresponding  stipulated  operational  values. 
Data  are  shown  in  both  logarithmic  and  physical  units.  It  should  be  remembered  that  the  critical 
region  of  the  profile  is  nearest  the  floor  of  the  test  section  in  the  vicinity  of  the  agent  and 
substrate.  Note  the  good  agreement  between  measured  and  operational  values. 

5.4  5-cm  Wind  Tunnel  Velocity  Profiles 

5.4.1.  Overview  of  Velocity  Profiles 


As  noted  previously,  a  design  goal  for  the  5-cm  Wind  Tunnel  was  to  allow  it  to  be 
replicated  in  large  numbers  (10-20)  and  used  in  multiple  locations.  Even  though  each  tunnel  had 
identical  dimensions,  it  was  important  to  demonstrate  that  each  of  the  tunnels  recreated  the 
desired  vertical  velocity  profiles.  Table  2  summarizes  the  extent  of  characterization  for  the 
initial  twelve  5-cm  Wind  Tunnels. 
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Figure  31.  Velocity  Profiles  Measured  at  the  Center  (Position  A)  of  the  5-cm  Wind  Tunnel 
in  Both  Physical  and  Log  Scales  (Tunnel  3B) 
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Table  2.  Summary  of  5-cm  Wind  Tunnel  Velocity  Profile  Characterizations 
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Position 

Tunnel 

Configuration 

Velocity 

A 

B 

C 

I> 

£ 

3A 

3 

Low 

X 

Medium 

X 

High 

3B 

3 

Low 

X 

X 

X 

X 

X 

Medium 

X 

X 

X 

X 

X 

High 

X 

X 

X 

X 

X 

3C 

3 

Low 

X 

Medium 

X 

High 

X 

3D 

3 

Low 

X 

Medium 

X 

High 

X 

3E 

3 

Low 

X 

Medium 

X 

High 

X 

3F 

3 

Low 

X 

X 

X 

X 

X 

Medium 

X 

X 

X 

X 

X 

High 

X 

X 

X 

X 

X 

3G 

3  Mod  1 

Low 

X 

Medium 

X 

High 

X 

3H 

3  Mod  1 

Low 

Medium 

High 
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3  Mod  1 

Low 

X 

Medium 

X 

High 

X 

3J 

3  Mod  1 

Low 

X 

Medium 

X 

High 

X 

3K 

3  Mod  1 

Low 

Medium 

X 

High 

3L 

3  Mod  1 

Low 

X 

X 

X 

X 

X 

Medium 

X 

X 

X 

X 

X 

High 

X 

X 

X 

X 

X 

The  original  goal  was  to  characterize  the  velocity  profiles  (i.e.,  measure  the 
velocity  profiles  for  all  three  test  conditions)  at  the  center  test  section  position  for  all  tunnels  and 
perform  a  complete  characterization  (at  all  five  test  section  positions  -  center  and  upstream, 
downstream,  left,  and  right  of  the  center)  for  at  least  one  each  of  the  Version  3  and  Version  3 
Mod  1  tunnels.  Two  repetitions  were  performed  for  each  test  condition.  At  the  time  of  this 
report,  complete  characterizations  have  been  completed  for  tunnels  3B,  3F  an  3L.  Tunnel  3D 
through  3F  are  expected  to  eventually  be  converted  to  the  Version  3  Mod  1  configuration  and 
will  be  re-characterized  at  the  center  location  only.  Tunnel  3H  has  not  been  characterized  as  yet. 

The  following  section  contains  the  plotted  vertical  velocity  profiles  for  all  of  the 
5-cm  Wind  Tunnels  along  with  the  corresponding  operational  profiles.  Profiles  are  shown  for 
the  positions  (Position  A  through  E)  and  velocities  tested  in  each  wind  tunnel  for  the  nominal 
35  °C  temperature.  The  results  of  additional  tests  to  evaluate  the  effect  of  temperature,  surface 
roughness  and  turbulence  screens  are  also  discussed.  While  vertical  velocities  were  measured 
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across  the  entire  height  of  the  test  section,  only  the  portion  of  the  profile  near  the  floor  is  of 
importance  since  this  is  the  region  that  affects  the  sorption  and  evaporation  rates  of  the 
agent/substrate  that  are  the  main  items  of  interest  in  the  Agent  Fate  Program.  Accordingly,  the 
turbulence  generators  were  intended  to  only  influence  the  bottom  portion  of  the  velocity  profile 
(i.e.,  between  2  cm  and  the  floor  of  the  test  section). 

Table  3  shows  the  difference  between  the  measured  and  operational  velocity 
values  at  a  height  of  1  and  2  cm  above  the  test  section  floor  for  all  of  tunnels  and  conditions 
tested  in  the  5-cm  Wind  Tunnel.  The  values  shown  are  the  average  of  the  repeated  runs.  These 
numbers  provide  an  indication  of  the  ability  of  the  5-cm  Wind  Tunnel  to  match  the  desired 
velocity  conditions  and  associated  velocity  profiles. 

At  2  cm,  the  velocity  values  at  the  center  of  the  piston  (position  A)  were  within 
about  3%  with  a  standard  deviation  of  about  12%  of  the  operational  value  for  all  of  the  tunnels. 
All  positions  had  slightly  higher  velocities  than  the  operational  value  ranging  from  3  to  12%. 

The  lateral  positions  (B  and  C)  had  differences  of  12  and  5%,  respectively,  with  about  the  same 
5%  deviation  with  the  position  closest  to  the  window  having  the  higher  values.  The  upstream 
and  downstream  positions  (D  and  E)  had  average  differences  of  about  4%  with  the  upstream 
position  having  the  largest  deviation  value  of  9%. 

The  values  at  I -cm  height  were  even  better.  The  velocities  at  the  center  of  the 
piston  (position  A)  were  within  about  1%  with  a  standard  deviation  of  about  4%  of  the 
operational  value  for  all  of  the  tunnels.  The  lateral  positions  (B  and  C)  were  similar  and 
exhibited  slightly  lower  velocities  ranging  from  2  to  6%  with  standard  deviations  of  about  4%. 
The  upstream  and  downstream  positions  (D  and  E)  had  average  differences  of  about  4  and  1%, 
respectively,  with  the  upstream  position  having  the  largest  deviation  value  of  16%. 

A  closer  examination  reveals  that  the  largest  average  values  and  standard 
deviations  are  those  at  the  2-cm  height  for  the  low  velocity  condition.  As  will  be  discussed  in 
Section  8,  this  represents  a  fully  laminar  profile  and  is  a  special  situation  as  compared  to  the  fully 
turbulent  profiles  of  the  medium  and  high  velocities.  For  the  latter  two  conditions,  the  average 
difference  and  standard  deviation  values  at  all  positions  and  both  heights  never  exceed  6%. 

These  observations  indicate  the  general  ability  of  the  5-cm  Wind  Tunnel  to 
accurately  and  repeatedly  achieve  the  reference  velocities  and  associated  velocity  profiles 
stipulated  by  the  Agent  Fate  test  matrix.  .  A  more  rigorous  quantitative  analysis  of  the  velocity 
profiles  based  on  boundary  layer  theory  is  contained  in  Section  8  of  this  report. 
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Table  3.  Percent  Difference  between  Measured  and  Operational  Velocities 
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5.4.2 


Tunnel  Version  3 


5.4.2. 1  Tunnel  3A 


Vertical  velocity  profiles  for  tunnel  3A  were  measured  only  at  the  low  and 
medium  velocities  at  the  center  position.  Characterization  was  cut  short  to  allow  the  tunnel  to  be 
installed  in  the  fume  hood,  instrumentation  incorporated  and  validation  tests  begun  so  that 
instrumentation,  testing  procedures  and  other  issues  could  be  identified.  As  can  be  seen  in 
Figures  32,  good  agreement  with  the  operational  profiles  was  achieved  with  very  good 
repeatability. 


Velocit>  (ni/s) 


Figure  32.  Measured  Velocity  Profiles  for  Tunnel  3A 


5. 4. 2. 2  Tunnel  3B 


While  tunnel  3A  was  being  used  to  evaluate  various  instrumentation  arrangements 
and  validate  the  5-cm  tunnel  design,  a  complete  velocity  characterization  was  performed  on 
tunnel  3B. 


Figure  33  shows  the  vertical  velocity  profiles  measured  at  the  center  position 
(Position  A)  for  all  three  velocities  compared  with  their  respective  operational  profiles  at  the 
center  position. 
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Figure  33.  Measured  Velocity  Profiles  for  Tunnel  3B  (Position  A) 


Figure  34  contains  the  vertical  velocity  profiles  measured  in  tunnel  3B  at  all  five 
positions  for  the  low  velocity  compared  with  the  operational  profile.  As  can  be  seen,  there  is 
virtually  no  difference  in  the  profile  across  the  piston  area. 
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Figure  34.  Tunnel  3B:  Low  Velocity 
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Figure  35  contains  the  vertical  velocity  profiles  measured  in  tunnel  3B  at  all  five 
positions  for  the  medium  velocity  compared  with  the  operational  profile.  The  only  notable 
deviation  from  the  operational  profile  occurs  at  location  B  (lateral  position  closest  to  the  test 
section  window)  where  the  velocities  are  slightly  lower  than  the  operational  profile.  However, 
even  at  this  position,  the  measured  profile  is  in  close  agreement  with  the  operational  profile  in 
the  critical  region  near  the  floor. 


Figure  35.  Tunnel  3B:  Medium  Velocity 


Figure  36  contains  the  vertical  velocity  profiles  measured  in  tunnel  3B  at  all  five 
positions  for  the  high  velocity  compared  with  the  operational  profile.  In  this  case,  the  largest 
deviation  from  the  operational  profile  occurs  at  positions  B  and  C  (lateral  locations)  where  the 
velocities  are  consistently  lower  than  the  operational  profile.  This  is  most  likely  due  to  side  wall 
effects.  However,  even  at  this  position,  the  measured  profile  is  in  close  agreement  with  the 
operational  profile  in  the  critical  region  near  the  floor.  At  a  height  of  1  cm,  the  difference  in 
velocity  is  on  the  order  of  5%  and  reduces  as  the  profile  approaches  the  test  section  floor. 
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Figure  36.  Tunnel  3B:  High  Velocity 


5. 4. 2. 3  Tunnel  3C 


Figure  37  contains  the  vertical  velocity  profiles  measured  in  tunnel  3C  at  the 
center  position  (position  A)  for  all  three  velocities  compared  with  their  corresponding 
operational  profiles.  Good  repeatability  is  present  as  well  as  close  agreement  with  their 
operational  profiles. 


— —  Low  Operational  ♦  Low  Measurement  1  A  Low  Measurement  2 

—  Medium  Operational  ♦  Medium  Measurement  1  A  Medium  Measurement  2 

- High  Theoretical  ♦  High  Measurement  1  A  High  Measurement  2 


Figure  37.  Measured  Velocity  Profiles  for  Tunnel  3C 
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5.4. 2.4  Tunnel  3D 


Figure  38  contains  the  vertical  velocity  profiles  measured  in  tunnel  3D  at  the 
center  position  (position  A)  for  all  three  velocities  compared  with  their  corresponding 
operational  profiles.  Good  repeatability  is  present  as  well  as  close  agreement  with  their 
operational  profiles. 


- Low  Operational  ♦  Low  Measurement  1  A  Low  Measurement  2 

——  Medium  Operational  ♦  Medium  Measurement  1  A  Medium  Measurement  2 

““  High  Operational  ♦  High  Measurement  1  A  High  Measurement  2 


Figure  38.  Measured  Velocity  Profiles  for  Tunnel  3D 


5.4. 2. 5  Tunnel  3E 


Figure  39  contains  the  vertical  velocity  profiles  measured  in  tunnel  3E  at  the 
center  position  (position  A)  for  all  three  velocities  compared  with  their  corresponding 
operational  profiles.  Good  repeatability  is  present  as  well  as  close  agreement  with  their 
operational  profiles. 
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Low  Operational  ♦  Low  Measurement  1  A  Low  Measurement  2 

—  Medium  Operational  ♦  Medium  Measurement  1  A  Medium  Measurement  2 

—  High  Operational  ♦  High  Measurement  1  A  High  Measurement  2 


Figure  39.  Measured  Velocity  Profiles  for  Tunnel  3E 


5. 4. 2. 6  Tunnel  3F 


Figures  40,  41,  and  42  contain  the  vertical  velocity  profiles  measured  in  tunnel  3F 
at  all  positions  for  the  low,  medium,  and  high  velocity,  respectively,  compared  with  their 
corresponding  operational  profiles. 
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-  Operational 
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♦  Position  H 
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■  Position  D 

■  Position  E 

Velocity  (m/s) 


Figure  40.  Tunnel  3F:  Low  Velocity 
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Figure41.  Tunnel  3F:  Medium  Velocity 


D  0.5  1  1.5  2  2.5  3  3.5  4 

Velocity  (m/i) 


Figure  42.  Tunnel  3F:  High  Velocity 
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5.4.3 


Version  3  Mod  1 


5.4.3. 1  Tunnel  3G 


Figure  43  contains  the  vertical  velocity  profiles  measured  in  tunnel  3G  at  the 
center  position  (position  A)  for  all  three  velocities  compared  with  their  corresponding 
operational  profiles.  This  particular  wind  tunnel  was  converted  from  a  Version  3  (horizontal 
transition  cone  and  turning  vane  section)  to  a  Version  3  Mod  1  (vertical  transition  cone  and 
turning  vane  section).  Velocity  profiles  were  measured  for  both  configurations.  While 
differences  in  the  measured  profiles  are  noted  between  the  two  tunnel  versions  at  the  larger  probe 
heights,  they  are  essentially  identical  in  the  critical  region  near  the  test  section  floor.  Good 
repeatability  is  present  as  well  as  close  agreement  with  their  operational  profiles. 


-  Operational 

♦  Tunnel  3G  -  Measurement  I 
A  Tunnel  3G  -  Measurement  2 

♦  Tunnel  3G,  Mod  I  -  Measurement  1 
■  Tunnel  3G,  Mod  I  -  Measurement  2 


Figure  43.  Measured  Velocity  Profiles  for  Tunnel  3G 


5. 4. 3. 3  Tunnel  31 


Figure  44  contains  the  vertical  velocity  profiles  measured  in  tunnel  31  at  the 
center  position  (position  A)  for  all  three  velocities  compared  with  their  corresponding 
operational  profiles.  Good  repeatability  is  present  as  well  as  close  agreement  with  their 
operational  profiles. 
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—  Operational 
♦  Measurement  1 
▲  Measurement  2 


Figure  44.  Measured  Velocity  Profiles  for  Tunnel  31 


5. 4. 3.4  Tunnel  3J 


Figure  45  contains  the  vertical  velocity  profiles  measured  in  tunnel  3J  at  the 
center  position  (position  A)  for  the  low  and  medium  velocities  compared  w  ith  their 
corresponding  operational  profiles.  Good  repeatability  is  present  as  well  as  close  agreement  w  ith 
their  operational  profiles. 
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—  Operational 
♦  Measurement  1 
▲  Measurement  2 


Figure  45.  Measured  Velocity  Profiles  for  Tunnel  3J 


5.43.5  Tunnel  3K 


Figure  46  contains  the  vertical  velocity  profiles  measured  in  tunnel  3K  at  the 
center  position  (position  A)  for  the  medium  velocity  compared  with  the  corresponding 
operational  profile.  As  noted  previously,  the  vertical  turning  vane  orientation  of  the  Version  3 
Mod  1  tunnel  required  the  addition  of  a  screen  at  the  entrance  to  the  fetch.  As  can  be  seen,  the 
screen  influences  the  profile  near  the  center  of  the  tunnel,  but  the  velocity  profile  near  the  test 
section  is  unaffected  and  shows  good  agreement  with  the  operational  profile. 
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Figure  46.  Tunnel  3K:  Medium  Velocity 
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5.4. 3. 6 


Tunnel  3L 


A  complete  velocity  profile  characterization  was  performed  for  tunnel  3L.  Figure 
47  contains  the  measured  velocity  profiles  in  tunnel  3L  for  all  five  positions  at  the  low  velocity 
condition.  As  was  the  case  with  the  Version  3  tunnels,  the  profiles  are  invariant  across  the  piston 
surface. 


0  0.5  1  1.5  2  2.5  3  3.5  4 

V  elocity  (m/s) 


Figure  47.  Tunnel  3L:  Full  Characterization  of  Low  Velocity 


Figure  48  contains  the  measured  velocity  profiles  for  all  five  positions  at  the 
medium  velocity  condition.  The  profiles  are  invariant  across  the  piston  surface  and  show  good 
agreement  with  the  operational  profile. 

Figure  49  contains  the  measured  velocity  profiles  for  all  five  positions  at  the  high 
velocity  condition.  Here,  the  profile  at  position  B  (lateral  position  furthest  from  test  section 
window)  has  slightly  lower  velocities  toward  the  center  of  the  tunnel  but  match  the  profile  near 
the  floor.  The  profile  for  Position  E  (downstream  of  the  piston  center)  near  the  floor  is  slightly 
different  from  the  operational  profile.  Again,  the  differences  in  velocities  are  only  on  the  order 
of  5%.  Otherwise,  the  profiles  are  invariant  across  the  piston  surface  and  show  good  agreement 
with  the  operational  profile. 
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Figure  48.  Tunnel  3L:  Full  Characterization  of  Medium  Velocity 
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Figure  49.  Tunnel  3L:  Full  Characterization  of  High  Velocity  Profile 
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5.4.4 


Effect  of  Temperature 


A  series  of  experiments  were  performed  in  Tunnel  3F  to  evaluate  the  effect  of  air 
temperature  on  the  vertical  velocity  profiles.  For  these  tests,  all  of  the  tunnel  heating  elements 
were  turned  off  which  essentially  produced  an  air  temperature  of  25  °C.  Measurements  were 
conducted  at  position  A  at  the  low  and  medium  velocities.  The  resulting  velocity  profiles  are 
shown  in  Figure  50  along  with  the  profiles  measured  for  the  nominal  35  °C  case  as  well  as  the 
operational  velocity  profiles.  As  can  be  seen,  the  air  temperature  had  minimal  effect  on  the 
profile. 


-  Operational 

■  Measured  -  35  deg.  C 
♦  Measured  -  25  deg.  C 


Figure  50.  Effect  of  Air  Temperature  on  Measured  Velocity  Profiles 

5.4.5  Effect  of  Surface  Roughness 

Tunnel  F  was  also  used  to  evaluate  the  effect  of  surface  roughness  on  the  resulting 
vertical  velocity  profile.  Two  different  concrete  substrates  were  evaluated,  each  measuring 

37.5  mm  in  diameter.  Both  had  the  same  surface  roughness  relatively  greater  than  glass.  The 
concrete  specimens  were  produced  with  a  surface  texture  depth,  to  meet  airfield  concrete 
specifications  for  hydroplaning  resistance.  Two  surface  texture  types  were  produced,  a  random 
brushed  texture  and  a  more  uniform,  ridged  texture.  The  concrete  studied  possessed  the  ridged 
surface  having  a  peak  to  trough  dimension  of  about  0.5  mm  and  a  peak  to  peak  dimension  of 
about  6.4  mm  resulting  in  about  6  ripples  across  the  diameter  of  the  substrate  samples.  Vertical 
profiles  were  measured  at  the  center  of  the  substrate  (position  A)  for  the  low  and  high  velocities. 
Tests  were  performed  with  the  concrete  wave  pattern  oriented  both  in-line  and  normal  to  the  air 
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flow  direction.  These  profiles  are  shown  in  Figure  5 1  compared  with  their  corresponding 
operational  profiles.  As  can  be  seen,  neither  the  increased  roughness  nor  the  orientation  of  the 
waves  had  any  appreciable  effect  on  the  profiles. 


““  Operational 
■  Measured  -Parallel  Grooves 

♦  Measured  -  Perpendicular  Grooves 


Figure  5 1 .  Effect  of  Concrete  Groove  Orientation  Relative  to  Flow  Direction 


5.4.6  Effect  of  Inlet  Screen 


Tunnel  K  was  the  first  Version  3,  Mod  1  tunnel  to  be  characterized.  As  noted 
previously,  the  vertical  orientation  of  turning  vanes  with  respect  to  the  vertical  velocity  profiles 
in  the  test  section  required  the  use  of  an  inlet  screen  at  the  entrance  to  the  fetch.  Tunnel  K  was 
used  to  develop  the  optimum  screen  configuration.  Various  screen  porosities  were  evaluated  and 
the  results  are  also  shown  in  Figure  52. 
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-  Operational  •  No  Screen  B  Window  Screen 

A  Medium  Screen  +  Ultra  Fine  Screen 


Figure  52.  Effect  of  Screen 

5.4.7  Consistency  between  Similar  Tunnel  Versions 

Figures  53  and  54  compare  vertical  velocity  profiles  measured  at  the  center 
position  in  five  different  Version  3  tunnels  and  three  different  Version  3  Mod  I  tunnels.  In  both 
cases,  good  consistency  is  achieved  for  both  tunnel  versions,  especially  in  the  critical  lower 
region  of  the  profiles  near  the  test  section  floor.  The  extreme  clustering  of  the  data  points  is 
indicative  of  the  good  repeatability  achieved  for  the  velocity  profiles  among  the  different  tunnels. 
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Probe  Height  (cm)  Probe  Height  (cm) 


—  Operational  Profile 

•  Measurement,  T unnel  3B  (Version  3) 

•  Measurement,  Tunnel  3C  (Version  3) 
A  Measurement,  Tunnel  3D  (Version  3) 

•  Measurement,  Tunnel  3E  (V  ersion  3) 

■  Measurement,  T unnel  3F  (V  ersion  3) 


Figure  53.  Consistency  of  Velocity  Profiles  in  Version  3  Tunnels 


-  Operational  Profile 

♦  Measurement,  Tunnel  3J  (V  ersion  3  Mod  1) 

A  Measurement,  Tunnel  3K  (V  ersion  3  Mod  1) 

•  Measurement,  Tunnel  3L  (Version  3  Mod  1) 


0  0.5  1  1.5  2  2.5  3  3.5  4 


Velocity  (m/s) 

Figure  54.  Consistency  of  Velocity  Profiles  in  Version  3  Mod  1  Tunnels 
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5.4.8 


Consistency  between  Different  Tunnel  Versions 


Figure  55  compares  the  velocity  profiles  measured  at  the  center  position  in  the 
Version  3  and  Version  3  Mod  I  configurations  of  the  5-cm  Wind  Tunnel.  As  can  be  seen,  good 
agreement  is  present  between  the  measured  and  operational  velocity  profiles  in  both  versions  of 
the  wind  tunnel. 


^  Operational 
♦  ■  Tunnel  3  (3F) 


A  •  Tunnel  3  Mod  1  (3L) 


Figure  55.  Velocity  Profiles  at  Center  Location  for  Tunnel  3  and  Tunnel  3  Mod  1 


6.  FREE  STREAM  TURBULENCE 

Free  stream  turbulence  is  the  random,  small  magnitude  velocity  fluctuations 
occurring  in  any  air  stream,  whether  in  an  atmospheric  boundary  layer  or  in  a  wind  tunnel.  The 
5-cm  Wind  Tunnel  forces  the  creation  of  a  turbulent  boundary  layer  by  introducing  turbulence  to 
the  flow  via  the  turbulence  generators  located  at  the  entrance  to  the  fetch  section  upstream  of  the 
test  section.  The  formation  of  a  turbulent  boundary  layer  and  its  characteristic  velocity  profile 
depend  on  having  a  certain  level  of  free  stream  turbulence  present  in  the  air  flow  (the  more 
turbulence,  the  fuller  the  profile).  However,  its  velocity  profile  (representing  a  specific  friction 
velocity  and  a  corresponding  shear  stress  at  the  surface)  is  relatively  insensitive  to  the  free 
stream  turbulence  level  as  long  as  it  remains  within  a  certain  range  of  values. 
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Free  stream  turbulence  intensities  were  measured  in  the  5-cm  Wind  Tunnel  for  all 
tests.  Representative  turbulence  intensities  for  each  of  the  three  velocity  conditions  are  shown  in 
Figure  56.  Flere  the  root  mean  square  of  the  three  orthogonal  velocity  fluctuations  (due  to 
turbulence)  is  referenced  to  the  local  total  velocity.  Thus,  while  the  absolute  free  stream 
turbulence  level  is  fairly  constant  throughout  the  tunnel,  the  intensity  increases  near  the  floor 
because  the  total  velocity  (the  denominator  in  the  turbulence  intensity)  decreases  near  the  floor. 
As  can  be  see,  the  turbulence  intensity  ranges  from  about  0.5  to  2%  for  the  low  velocity  and 
about  4  to  1 4%  for  the  medium  and  high  velocities. 

The  turbulence  level  for  the  low  velocity  profiles  indicates  a  laminar  velocity 
profile.  The  flow  lacks  the  inertia  energy  to  overcome  the  viscous  effects  and  develop  a  fully 
turbulent  boundary  layer.  On  the  other  hand,  the  medium  and  high  velocity  profiles  contain 
sufficient  inertia  energy  to  develop  a  fully  turbulent  boundary  layer  profile.  However,  the  free 
stream  turbulence  level  in  the  wind  tunnel  does  not  affect  the  evaporation  rates  of  the  agents  as 
long  as  the  friction  velocities  are  compatible." 


Figure  56.  Turbulence  Intensity  of  5-cm  Wind  Tunnel 
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7. 


EXPERIMENTAL  OPERATIONS 


The  test  matrix*  for  the  5-cm  Wind  Tunnel  includes  the  agents  and  substrates  as 
well  as  the  test  parameters  consisting  of  temperature,  drop  size,  velocity  and  relative  humidity. 

The  agents  to  be  investigated  are  shown  listed  in  order  of  their  priority,  based  on 
importance  to  the  program  and  feasibility  of  the  measurements  given  their  complexity: 

HD 

VX 

GD 

Thickened  VX 

Thickened  GD 

A  key  feature  of  the  substrate  was  its  interaction  with  the  agent.  The  substrates 
were  resolved  into  four  basic  categories: 

1 .  non-absorptive,  non-reactive 

2.  absorptive,  non-reactive 

3.  non-absorptive,  reactive 

4.  absorptive,  reactive 

Each  agent  might  have  specific  substrates  assigned  to  each  category,  depending 
mostly  on  its  reactivity.  The  specific  substrates  to  be  tested  were  prioritized  as  follows  along 
with  their  general  substrate  interaction  characteristics  noted  above: 


Substrate 


Cateuory 


1.  Glass 

2.  Dry  Sand/Clay 

3.  Concrete 

4.  Asphalt 

5.  Grass 


non-absorptive,  non-reactive 
absorptive,  non-reactive 
absorptive,  reactive 
absorptive,  reactive 
absorptive,  reactive 


A  series  of  validation  tests  were  performed  involving  HD  on  glass  (glass 
representing  a  non-absorbent,  non-reactive  substrate  for  HD).  With  the  “Validation"  tests 
completed,  the  program  then  entered  the  “production"  testing  phase.  This  involves  testing  agents 
with  the  more  complex  substrates.  The  initial  production  tests  involved  HD  on  sand.  Each 
combination  of  agent  and  substrate  were  tested  at  each  of  the  three  velocity  conditions,  three 
temperatures  and,  if  relevant,  three  relative  humidity  values  as  presented  in  Figure  57. 


’  Kilpatrick,  W.  T.;  Ling,  E.  E.;  Hin.  A.  R.  T.;  Brevett.  C.  A.  S.;  Fagan,  M.  W.;  Murdock.  W.  P.  Jr.  Testing 
Requirements  for  Predictive  Model  Development  Using  Environmental  Wind  Tunnels.  Air  Force  Research 
Laboratory:  Wright-Palterson  Air  Force  Base.  OFI.  submitted  for  publication. 
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Parameter  Code  Kev: 


Drop  Size 
(MU 


Figure  57.  Basic  Test  Matrix  Arrangement  and  Coding* 


8.  EVALUATION  OF  5-CM  WIND  TUNNEL  VELOCITY  PROFILES 

8.1  Overview  Evaluation 


An  analysis  was  performed4  to  evaluate  how  well  the  experimental  velocity 
profile  measurements  in  the  5-cm  Wind  Tunnel  agree  with  the  operational  profiles  specified  by 
the  Agent  Fate  program  and  the  effect  of  any  differences  on  the  subsequent  agent  evaporation 
measurements. 


8.2  Turbulent  Boundary  Laver,  Log-Law,  and  Law-of-the-Wall  Variables 

Fully  developed  boundary  layers  or  internal  flows  include  a  region  identified  as 
the  logarithmic-layer  where  the  velocity  profile  can  be  represented  by: 


u 

Ut 


u  v 

=  — ln(— ^-)  +  C 

K  V 


(1) 


Kilpatrick,  W.  T.;  Ling,  E.  E.;  Hin,  A.  R.  T.;  Brevelt,  C.  A.  S.;  Fagan,  M.  W.;  Murdock.  W.  P.  Jr.  Testing 
Requirements  for  Predictive  Model  Development  Using  Environmental  Wind  Tunnels.  Air  Force  Research 
Laboratory:  Wright-Palterson  Air  Force  Base,  OH,  submitted  for  publication. 
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Where: 


u  =  velocity  parallel  to  the  surface 
y  =  distance  normal  to  the  surface 
v  =  kinematic  viscosity 

k  =  von  Karman  Universal  Constant  of  Turbulence  (  k  =  0.4 ) 


uT  =  friction  velocity  (  u  t 


Tw  =  stress  acting  on  surface 
p  =  air  density 

C  =  5.5  (fully  developed,  turbulent  boundary  layer  on  a  flat  plate  with  a  zero 
pressure  gradient) 


A  Fortran  program  was  developed  to  take  three  or  more  data  points  from  the 
measured  profile  and  to  perform  a  linear  least  squares  fit  to  u  versus  ln(y)  to  compute  the  slope 
and,  assuming  k  =  0.4,  to  determine  the  most  likely  value  of  ut  and  the  intercept  C.  The  initial 

points  were  selected  by  trial  and  error  but  usually  excluded  the  2  or  3  points  nearest  the  floor  and 
the  points  beyond  the  centerline  of  the  tunnel  which  deviated  from  the  linear.  Most  of  the 
primary  profile  measurements  were  made  at  the  center  of  the  test  section  (Position  A)  and  were 
repeated  up  to  three  times  and  these  runs  were  grouped  together.  The  numerical  data  were 
transferred  to  a  plotting  routine  where  the  least  square  fits  were  compared  with  the  data. 


As  noted  previously,  a  set  of  operational  velocity  profiles  were  defined  for  the 
Agent  Fate  Wind  Tunnels  representing  wind  velocities  of  0.5,  3.0,  6.0  m/s  at  a  height  of  2  m 
above  the  ground.  It  is  assumed  that  the  smooth  wall  equilibrium  turbulent  boundary  layer  exists 
within  the  2  m.  The  previous  equation  was  used  to  determine  the  associated  friction  velocities 
and  the  velocities  at  a  2-cm  height  for  the  three  operational  velocity  profiles  with  k  =  0.4  and 
C  =  5.5. 


Table  4.  Extrapolated  Operational  Velocities  to  2-cm  Height 


Velocity  (m/s)  @  2  m 

ut.  Friction  Velocity  (m/s) 

Velocity  (m/s)  @  2  cm 

0.5 

0.020 

0.26 

3.0 

0.1038 

1.79 

6.0 

0.1966 

3.70 
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The  boundary  layer  profile  can  also  be  expressed  in  terms  of  no-dimensional 
“law  of  the  wall”  coordinates  (u  and  y+): 


Where: 


u+  =  Law  of  the  wall  velocity  coordinate  =  u/uT 

uT  =  Friction  or  shear  velocity  =  Vu(<9 u/<9y)w 
y+  =  Law  of  the  wall  height  coordinate  =  uTy/u 

(du/dy)w  =  velocity  gradient  at  surface  (y  =0) 

A  hot  wire  anemometer  was  used  to  measure  the  velocity  profiles  and  the 


turbulence  intensity  distribution  in  the  5-cm  Wind  Tunnels.  These  include  two  slightly  different 
configurations:  one  having  a  horizontal  inlet  and  the  other  having  a  vertical  inlet.  Velocity 
profiles  were  measured  at  the  center  of  the  test  section  for  three  velocity  conditions.  These  are 
designated  as  high,  mid  and  low  for  the  2  m  height  values  for  the  operational  velocities  of  6,  3, 
and  0.5  m/s,  respectively.  Selected  points  are  shown  in  “law  of  the  wall”  coordinates  in 
Figure  58.  As  can  be  seen,  for  the  low  speed  case,  the  range  of  the  measurements  are  within  the 
transition  region  and  do  not  extend  into  the  log-law  region.  Whereas  the  high-speed  case  is 
entirely  within  the  log-law  and  the  mid  speed  is  predicted  to  have  a  substantial  range  in  the  log- 
law  but  with  the  lower  points  extending  into  the  transition  region. 
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Figure  58.  Non-dimensional  Velocity  Profile 
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8.3 


Mid  and  High  Velocities 


Two  series  of  profile  measurements  are  considered  representing  the  horizontal 
inlet  configuration  with  no  inlet  screen  (T3B)  and  the  vertical  inlet  configuration  with  the 
standard  44%  porosity  screen  (T3L).  The  friction  velocity  defines  the  near-wall  velocity 
distribution  for  a  smooth  wall  and  the  specification  of  the  velocity  at  2  m  along  with  the  assumed 
profile  form  (Log-Law)  defines  the  friction  velocity.  The  friction  velocity  was  determined  from 
the  slope  of  u  versus  ln(y)  assuming  the  von  Karman  universal  constant  of  turbulence  k  =  0.4. 

A  linear  least  squares  technique  was  used  with  3-7  data  points  identified  as  falling  in  the  linear 
range.  Figure  59  illustrates  the  quality  of  the  fit. 

Figures  60  and  61  present  the  measured  velocity  profile  values  in  non- 
dimensional  form  compared  with  the  operational  profiles.  As  can  be  seen,  in  both  cases,  the 
measurements  are  higher  than  the  operational  values  but  the  slopes  are  in  general  agreement. 
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Figure  59.  Log-Law  Slope  (Tunnel  T3L) 
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Figure  60.  Measured  vs.  Operational  Velocity  Profile  Values  of  Tunnel  T3B 
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Figure  61.  Measured  vs.  Operational  Velocity  Profile  Values  of  Tunnel  T3L 


0.5  m/s 

A 


6  m/s 

3  m/s  ^ 


-  Sublayer  Transition 
Log-Law 

A 

Ref.  Pts 

■ 

Range  of  Exp  High. 

□ 

Range  of  Exp.  Mid 

H 

Range  of  Exp.  Low 

-  IT=  Y 

o 

T3L-MA1-2 

♦ 

T3L-HA1-2 

The  shift  in  the  measurements  relative  to  the  operational  curve  is  reflected  in  the 
intercept  C  obtained  from  the  least  square  fit.  Since  there  is  no  trend  with  velocity  in  these  data, 
an  overall  average  of  C  =  8.5  ±  1 .0  describes  the  data  whereas  the  assumed  value  for  the 
operational  profile  is  C  =  5.5.  There  could  be  several  causes  for  this  modification  of  the  profile 
from  the  normal  including:  a  departure  from  the  two-dimensionality  condition  because  of 
sidewall  effects  due  to  the  small  scale  of  the  facility;  the  inlet  flow  moving  through  a  90°  bend 
prior  to  entering  the  fetch;  and  the  unsymmetrical  surface  roughness  created  by  the  turbulence 
generators  located  only  five  test  section  heights  upstream  of  the  test  section. 
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Not  all  the  profiles  investigated  show  such  a  large  value  for  the  intercept  C.  A 
more  complete  picture  can  be  gained  from  Table  1  where  friction  velocity  and  C  are  presented 
for  fifteen  profiles  measured  in  different  5-cni  Wind  Tunnels.  The  average  of  all  the  runs  is 
C  =  6.6  and  u  x  =  0. 1 80  and  0.096  for  the  high  and  mid  velocities  respectively.  Using  these 
results  it  is  possible  to  work  back  to  obtain  the  velocity  at  2  m  for  the  average  measured  profile, 
i.e.,  u,m  =  5.7  and  2.9  m/s,  which  differ  only  slightly  from  the  operational  values. 

Table  5.  Friction  Velocity  Data  Obtained  from  Measured  5-cm  Wind  Tunnel  Profiles 


RUNS 

11 1  m/s 

2  a  ut  % 

c 

COMMENT 

T3B-HA1-3 

0.160 

23 

9.2 

Horiz.  Inlet,  No  Screens 

T3F-HA9-10 

0.202 

6.0 

4.5 

Horiz.  Inlet,  No  Screens 

T3G-HA3-4 

0.189 

3.7 

5.2 

Horiz.  Inlet,  No  Screens 

T3F-HAI 1-12 

0.188 

6.0 

5.8 

Horiz.  Inlet,  No  Screens 

AVERAGE 

0.184 

6.2 

T3L-HA1-2 

0.174 

8.3 

7.2 

Vert.  Inlet,  44%  Porosity  Screens 

Table  5.  Friction  Velocity  Data  Obtained 
Tunnel  Profiles  (Conti 

rom  Measured  5-cm  Wind 
inued) 

RUNS 

ur  m/s 

2  a  u  t  % 

C 

COMMENT 

T3EYE-HA3-4 

0.201 

4.3 

4.7 

Vert.  Inlet,  44%  Porosity  Screens 

T3K-HA8-9 

0.158 

6.4 

9.4 

Vert.  Inlet,  44%  Porosity  Screens 

AVERAGE 

0.178 

6.9 

T3G-MA1-2 

0.096 

5.4 

6.1 

Horiz.  Inlet,  No  Screens 

T3F-MA5-6 

0.103 

9.3 

5.2 

Horiz.  Inlet,  No  Screens 

T3F-MA7-8 

0.108 

7.6 

4.3 

Horiz. Inlet,  No  Screens 

T3B-MA 1-3 

0.088 

17 

7.9 

Horiz.  Inlet,  No  Screens 

T3F2-MA2 

0.102 

13.1 

5.4 

Horiz.  Inlet,  No  Screens 

AVERAGE 

0.100 

5.6 

T3K-MA6-7 

0.099 

6.8 

5.8 

Vert.  Inlet,  44%  Porosity  Screens 

T3L-MA1-2 

0.078 

8.3 

10.2 

Vert.  Inlet,  44%  Porosity  Screens 

T3EYE-MAI-2 

0.099 

5.6 

6.08 

Vert.  Inlet,  44%  Porosity  Screens 

AVERAGE 

0.092 

7.4 
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Figure  62  presents  the  T3B  experimental  data  in  physical  variables  and  a  semi- 
logarithmic  format  and  includes  the  operational  profile  which  is  computed  based  on  the  mixing 
length  analysis  the  u 4  =  y +  line  referred  to  the  friction  velocity. 
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Figure  62.  Semi-Logarithmic  T3B  Profiles  vs.  Operational  Profile 


The  velocity  at  a  height  of  2  cm  was  initially  used  as  the  reference  value  for  the 
5-cm  Wind  Tunnel  validation  tests  involving  HD  on  glass.  The  reference  velocity  for  the  follow- 
on  production  tests  was  later  changed  to  that  at  a  1-cm  height.  Table  5  shows  how  well  the 
velocities  measured  at  both  reference  heights  in  four  different  5-cm  Wind  Tunnels  compared 
with  the  associated  operational  values. 

Thus,  the  measured  velocity  at  a  2-cm  height  is  about  5-7%  lower  than  the 
operational  velocity  with  an  average  difference  of  3.9%  and  a  standard  deviation  of  1 .3%.  The 
measured  velocity  at  an  1-cm  height  4  to  5%  lower  than  the  operational  velocity  with  an  average 
of  4.0%  with  a  standard  deviation  of  1 .3%.  The  least  square  fit  u,cm  agrees  closely  with  the 
experimental  data  which  again  shows  that  the  data  and  fit  diverge  somewhat  from  the  operational 
profile. 
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Table  6.  Comparison  of  Measured  and  Operational  Velocities  at  1-  and  2-cm  Heights 


Tunnel/Run 

Operational 
Profile  (m/s) 

Measured 
Profile  (m/s) 

Difference 

% 

LS  Fit 

Profile  (m/s) 

Velocities  at  1  cm 

T3B-HA 1-3 

3.45 

3.37 

-2.2 

3.32 

T3L-HA7-8 

3.45 

3.33 

-3.6 

3.30 

T3B-MA1-3 

1.66 

1.58 

-4.7 

1.58 

T3L-MA1-2 

1.66 

1.57 

-5.2 

1.59 

Velocities  at  2  cm 

T3B-HA1-3 

3.69 

3.66 

-0.8 

3.57 

T3L-HA7-8 

3.70 

3.57 

-3.6 

3.62 

T3B-MA1-3 

1.80 

1.71 

-5.0 

1.74 

T3L-MA1-2 

1.79 

1.67 

-6.7 

1.71 

8.4  Low  Speed  Velocity 

The  low  velocity  data  would  be  expected  to  fall,  for  the  most  part,  in  the  transition 
or  buffer  region  between  the  laminar  sub  layer  and  the  fully  turbulent  or  the  log-law  region. 

Thus,  the  analysis  of  the  [low  velocity]  profile  assuming  a  semi-logarithmic  region  is  not 
applicable  as  it  is  for  the  higher  [two]  velocities.  The  measurements  extended  sufficiently  near 
the  wall  so  that  the  velocity  gradient  at  the  wall  is  obtained  by  extrapolation  and  thus  the  wall 
shear  stress  and  the  friction  velocity  can  be  determined.  The  existing  computer  code  was 
modified  to  skip  the  log-law  analysis  and  perform  the  extrapolation  in  the  following  form 

u  =  Ay  +  By2  +  Cy3  (2) 

In  this  form  the  boundary  condition  u  — >  0  as  y  — >  0  is  automatically  satisfied. 
To  actually  carryout  the  fitting  of  the  data  to  determine  the  coefficients,  the  data  were  put  into 
the  form  of  u  /  y  vs.  y  or 

—  =  A  +  By  +  Cy2  (3) 

y 

The  results  showed  that  the  data  point  closest  to  the  wall  had  to  be  excluded  in 
order  to  obtain  anything  like  a  monotonic  fit  and  then  only  a  limited  number  of  profile  data 
points  could  be  used.  The  shape  of  the  measured  profile  is  more  complicated  than  can  be 
described  by  such  a  simple  formula  without  introducing  inflections  in  the  fitted  curve.  The 
resulting  T3L  and  T3B  data  are  shown  in  Table  6. 
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Table  7.  Low  Velocity  Profile  Analysis  Results 


Tunnel 

From  Data 

Operational 

ut  (m/s) 

8(m) 

us  (m/s) 

ut  (m/s) 

T3L 

0.0396 

0.011 

0.221 

0.02 

T3B 

0.0390 

0.010 

0.216 

0.02 

The  ability  of  the  quadratic  fit  to  describe  the  near  wall  points  for  the  T3L  case  is 
shown  in  Figure  63  compared  with  the  operational  profile. 


Figure  63.  Experimental  vs.  Operational  Low  Speed  Velocity  Profiles 


The  friction  velocity  calculated  by  the  above  procedure  is  a  factor  of  two  larger 
than  the  operational  profile.  It  is  difficult  to  see  how  the  existing  data  could  be  extrapolated  to 
the  wall  with  a  wall  slope  as  small  as  required  (by  a  factor  of  2.)  Linearly  extrapolating  the 
next  to  the  last  data  point  gives  a  value  of  ut  =  0.03  that  is  still  50%  larger  than  the  Operational 
prediction.  It  is  possible  that  the  hot  wire  anemometer  data  at  this  low  velocity  could  be 
contaminated  by  wall  interference.  Assuming  the  uT  obtained  in  this  way  allows  the  profile  to 
be  plotted  in  the  law  of  the  wall  coordinates  as  shown  Figure  64.  In  these  coordinates  the  data 
fall  well  below  the  sub  layer  u  '  =  y  line  and  the  empirical  buffer  layer  line. 
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18 


Figure  64.  Low  Velocity  “Law  of  the  Wall”  Coordinates 


The  deviation  from  the  buffer  layer  distribution  in  the  law  of  the  wall  coordinates 
does  not  correspond  to  the  desired  turbulent  profile.  It  appears  appropriate  to  consider  the  low 
speed  data  to  be  a  laminar  boundary  layer  with  some  imposed  turbulence  from  the  tripping 
elements.  The  friction  velocity  correlates  with  that  of  a  laminar  boundary  layer  just  as  well  as 
for  the  turbulent  boundary  layer.  Thus,  u  t  represents  the  characteristic  convection  velocity  for 

the  laminar  evaporation  just  as  well  as  for  the  turbulent  flow.  For  the  same  ut  a  laminar 

boundary  layer  may  be  expected  to  produce  the  same  evaporation  rate  as  long  as  the 
concentration  distribution  in  the  vicinity  of  the  drop  is  not  affected  by  the  turbulence  away  from 
the  fluid  surface.  Thus,  the  low  speed  velocity  condition  is  quite  different  from  that  of  the  higher 
two  velocities  and  may  be  considered  to  represent  a  “very  low”  wind  speed  case. 
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8.5 


Turbulence  Intensity  Measurements 


Turbulent  intensity  measurements  were  obtained  in  the  form  of  root-mean-square 
stream-wise  fluctuations.  In  order  to  emphasize  the  sub  layer,  the  data  are  re-normalized  by  the 
friction  velocity  rather  than  the  local  velocity.  The  turbulent  intensity  data  were  converted  using 

u’+  =  —  u+  and  the  results  for  selected  profiles  are  compared  in  Figure  65.  The  T3B  data  refer  to 
u 

a  horizontal  inlet  tunnel  without  screens  while  the  vertical  inlet,  T3K  tunnel  has  inlet  screens  of 
61  and  44%  porosity.  As  might  be  expected  the  turbulence  intensity  decreases  with  decreasing 
screen  porosity. 
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Figure  65.  Comparison  of  Turbulence  Intensities 


Figure  66  provides  examples  of  the  wall  layer  turbulence  level  for  all  three 
velocities.  The  maximum  of  the  low  velocity  intensity  data  normalized  by  the  friction  velocity  is 
more  than  a  factor  of  100  less  than  the  mid  and  high  velocity  cases  for  most  of  the  runs.  This 
provides  additional  evidence  that  the  low  velocity  data  should  be  considered  laminar  and  treated 
differently  than  the  higher  velocity  data. 


73 


3.5  I 

3 

2.5 

2 

1.5 

1 

D.5 


A 

A 


— • — 

KlebanofT 

A 

Klewicki  A  Falco 

■ 

fi  wire  pm  bp 

♦ 

T3B  HA1  3 

■ 

T3L  HA7 

A 

T3B  MAI 

A 

T3L  MAI 

♦ 

T3B  LAI 

□ 

T3K-LA1  2 

A 

T3L  LA13 

Estimate 

-E - £l— 

1D 


♦  ♦  ♦  ♦  ♦  ♦  ♦♦♦, 


□□DO 


1DD 


10Q0 


Y* 

Figure  66.  Effect  of  Velocity  on  Turbulence  Intensity 


8.6  Effect  of  Differences  between  Measured  and  Operational  Velocity  Profiles 

The  analysis  provides  a  quantitative  comparison  of  how  well  the  velocity  profiles 
measured  in  the  5-cm  Wind  Tunnel  agree  with  the  specified  operational  values.  At  the  high  and 
mid  velocity  conditions,  corresponding  to  the  operational  profiles  at  6  and  3  m/s,  respectively, 
the  friction  velocities  are  0.180,  0.096  m/s  with  an  average  intercept  C  of  6.6  which  are  different 
from  the  Operational  values  of  0. 1966,  0. 1038  m/s  and  C  of  5.5.  Thus,  the  experimental  friction 
velocity  data  are  consistently  lower  by  7.5  and  8.4%,  respectively  than  specified  by  the  reference 
and  the  intercept  is  consistently  higher  by  1 8%. 

A  more  meaningful  criterion  to  quantify  this  comparison  is  to  relate  it  to  its 
subsequent  effect  on  the  agent  evaporation  rate  which  represents  the  ultimate  purpose  of  the 
5-cm  Wind  Tunnel.  A  two-dimensional  model  was  formulated  based  on  a  Couette  flow  that 
showed  that  the  characteristic  convection  velocity  is  the  friction  velocity.  The  effect  of 
differences  between  the  experimental  and  operational  velocity  profiles  can  be  compared  in  terms 
of  how  well  their  respective  friction  velocities  agree.  The  slope  of  the  turbulent  semi-logarithmic 
profile  is  used  to  determine  the  experimental  friction  velocities  from  the  measurements.  The 
uncertainties  in  individual  tests  in  terms  of  two  standard  deviations  in  the  friction  velocity  vary 
considerably  with  a  high  of  23%  and  a  low  of  4%  with  an  overall  average  of  8.7%.  Thus  the 
potential  uncertainty  is  of  the  same  order  as  the  deficit  in  the  friction  velocity  with  respect  to  the 
operational  profile.  The  evaporation  model  modified  by  initial  test  data  predicts  the  evaporation 
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rate  to  depend  on  the  friction  velocity  to  the  1/2  power  and  thus  the  uncertainty  in  rate  due  to 
friction  velocity  could  be  of  the  order  of  4-5%,  on  average. 

The  intercept  in  the  Log-Law  varies  considerably  in  individual  tests  from  a  high 
of  just  over  10  to  a  low  in  the  4’s  with  an  average  of  6.6.  In  principle,  the  intercept  is  not 
relevant  to  the  sub  layer  velocity  gradient  and  friction  velocity.  However  it  is  of  concern  in  that 
it  is  an  indication  of  the  extent  of  the  sub  layer  and  the  rapidity  of  the  damping  of  the  turbulence 
as  the  wall  is  approached.  In  addition,  it  may  be  an  indication  of  the  lack  of  two-dimensionality 
of  the  flow. 


The  low  velocity  condition  is  a  special  case.  Theoretically  the  operational  profile 
predicts  that  most  of  the  low  speed  data  to  be  mostly  in  the  buffer  layer  and  that  the  semi- 
logarithmic  slope  cannot  be  used  to  determine  the  fiction  velocity.  As  an  alternative,  the  data 
were  extrapolated  to  the  wall  and  the  slope  of  the  resulting  curve  used  to  determine  the  shear 
stress  and  friction  velocity.  The  results  show  the  friction  velocity  to  be  0.04  m/s  compared  to  the 
expected  value  of  0.02  m/s.  The  data  plotted  in  law  of  the  wall  coordinates  show  that  the 
experimental  data  deviate  significantly  from  that  of  a  turbulent  profile,  suggesting  the  low  speed 
measurements  were  essential  laminar.  Examination  of  the  turbulent  intensity  data  confirm  that 
the  velocity  fluctuation  level  is  substantially  lower  than  observed  for  the  other  velocities.  It  is 
probably  impossible  to  produce  a  turbulent  boundary  layer  at  such  a  low  velocity  in  a  5-cm 
tunnel.  The  friction  velocity  of  the  laminar  flow  can  also  be  related  to  a  hypothetical  2  m 
turbulent  boundary  layer  at  a  velocity  of  about  1  m/s. 


9.  CONCLUSIONS 

1 .  A  special  5-cm  Wind  Tunnel  was  designed  specifically  to  measure  the  release 
and  retention  of  Chemical  Warfare  Agents  (CWAs)  from  various  materials  under  constant, 
simulated  environmental  conditions.  These  wind  tunnels  and  associated  instrumentation  provide 
a  new  capability  for  measurement  of  the  environmental  fate  of  any  toxic  chemical  contaminant 
on  most  surfaces. 

2.  A  wind  tunnel  was  successfully  designed  that  meets  the  requirements  for  a 
reproducible  and  uniform  flow  field  around  contaminant  drops  on  a  surface  and  the  tunnel  can  be 
fitted  within  the  dimensions  of  a  standard  chemical  fume  hood  that  meets  safety  requirements  for 
CWA  operations. 

3.  The  tunnel  can  create  vertical  velocity  profiles  that  match  the  lower  portion 
of  a  wind  induced  atmospheric  boundary  layer. 

4.  The  specified  velocity  profiles  are  uniform  across  the  area  containing  the 
droplet(s)  on  surface  materials,  within  the  wind  tunnel  test  section. 

5.  The  specified  velocity  profiles  can  be  characterized  by  a  single  reference 
velocity  at  a  specified  height  in  the  test  section  and  this  reference  velocity  can  then  be  employed 
during  each  experiment  to  set  and  then  continuously  monitor  the  wind  speed. 
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6.  Free  stream  turbulence  in  the  5-cm  Wind  Tunnel  has  similar  values  to  those 
of  a  typical  atmospheric  boundary  layer. 

7.  The  wind  speed  profile  measurements  and  fit  to  the  operational  profile  are 
closely  matched  and  reproducible;  therefore,  the  hot  wire  measurements  during  a  wind  tunnel 
experiment  can  be  related  to  a  wind  speed  measured  at  a  2-m  height  or  other  height  at  a  weather 
station  at  an  airfield  or  other  Fixed  Site.  The  2-m  height  wind  speed  data  can  then  be 
continuously  input  into  a  real  time  model  during  a  chemical  attack  and  the  volatilization  rate  can 
be  estimated  from  the  5-cm  Wind  Tunnel  data  in  which  the  wind  speed  was  monitored  at 

10  mm,  or  any  similar  reference  height. 
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APPENDIX  I  -  HOT  WIRE  ANEMOMETER  CALIBRATION  PROCEDURE 


The  calibration  of  a  single  component  hot  wire  anemometer  is  relatively  straight 
forward.  Typically  better  sensitivity  can  be  achieved  if  the  calibration  is  limited  to  the  range  of 
velocities  expected  in  the  experiments.  For  the  5-cm  Wind  Tunnel  this  wind  speed  range  was 
selected  to  be  approximately  0.05  to  10-m/s.  Knowing  the  calibration  wind  speed  range,  the 
atmospheric  pressure  and  room  air  temperature,  a  calibration  spreadsheet  was  generated  that 
provided  the  different  tunnel  configurations  and  pressures  required  to  achieve  each  wind  speed. 
The  calibration  spread  sheet  is  shown  in  Table  ALL 

Table  AI-1.  Hot  Wire  Anemometer  Calibration  Spreadsheet 


Atm  Press 
Atm  Temp 
K 


757  mm-hg  Enter  Atm  Press  and  Temp  on  "No  Plate  Sheet" 

26.1  deg  C 

1 .021023  Standard  Conditions  29.92  in-Hg  (760  mm-Hg)  and  70  deg  F 


Vetoed 

y 

(mfr) 

Plate  82 

Plate  81 

No  Plate 

(in 

H20) 

0.1"  WIG 
(NFS) 

0.5" 

WIG 

(NFS) 

5.0" 

WIG 

(NFS) 

(in 

H20) 

5.0" 

WIG 

(NFS) 

(in 

H20) 

0.5- 

WIG 

(NFS) 

5.0"  WIG 
(NFS) 

1 

0.05 

0.0056 

5.6 

2 

0.01 

0.0130 

13.0 

3 

0.11 

0.0301 

30.1 

4 

0.17 

0.0694 

69.4 

5 

0.25 

0.1602 

32.0 

6 

0.38 

0.3698 

74.0 

7 

0.57 

0.8536 

17.1 

1 

0.15 

1.9705 

39.4 

9 

1.28 

4.5490 

91.0 

10 

1.92 

0.263 

5.3 

11 

2.88 

0.597 

11.9 

12 

4.32 

1.355 

27.1 

13 

6.49 

3.079 

61.6 

14 

9.73 

0.223 

44.6 

15 

14.60 

0.500 

10.0 

16 

21.89 

1.121 

22.4 

17 

32.84 

2.514 

50.3 

IS 

45.00 

4.708 

94.2 

Using  the  atmospheric  pressure  and  room  temperature  as  input,  a  correction  factor 
K  is  calculated  to  determine  the  true  air  density  for  non-standard  conditions.  The  spreadsheet  is 
used  to  compute  the  difference  between  the  test  section  and  settling  chamber  pressures  which  is 
then  used,  in  combination  with  the  appropriate  orifice  plate  and  electronic  pressure  gage 
required,  to  achieve  the  desired  calibration  velocity.  Once  these  setting  have  been  established, 
the  calibration  can  be  performed. 

Another  important  aspect  of  hot  wire  anemometry  that  should  be  noted  is  the 
effect  of  different  air  temperatures  on  the  resulting  measurements.  Because  hot  wire/film 
anemometry  relies  on  the  principles  of  heat  transfer,  a  difference  between  the  calibration  air 
temperature  and  the  air  temperature  of  the  actual  experiments  can  introduce  an  error  of 
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approximately  1%  per  °C  of  difference.  The  test  matrix  temperatures  for  the  Agent  Fate  program 
consisted  of  three  constant  values  ( 15,  35  and  50  °C).  Because  of  the  calibration  instrumentation 
available,  a  room  temperature  of  23-25  °C  was  the  typical  calibration  temperature.  Thus,  an 
ambient  air  temperature  difference  occurred  between  the  calibration  and  velocity  profile 
measurements  which  required  a  correction  to  be  performed  on  the  resulting  measurements. 

An  analytical  correction  is  the  most  popular  and  widely  used  temperature 
compensation  technique  and  relies  on  the  principles  of  heat  transfer  and  requires  a  combined 
measurement  of  both  the  fluid's  velocity  (i.e.,  voltage  difference)  and  temperature.  Before  the 
output  voltage  is  converted  to  a  velocity,  the  following  expression  is  applied  to  correct  for  the 
temperature  difference. 


Ecorrcclcd  E0ui  [(Ts  T|Cst)/(Ts  Tca|)] 


12 


Where: 

Ecorrcctcd  =  corrected  bridge  output 

Ts  =  hot  wire/film  sensor  temperature  (typically  250  °C) 

Tca|  =  calibration  temperature  (23  to  25  °C  for  this  study) 

T.es,  =  air  temperature  during  the  test 

Eou,  =  uncorrected  bridge  output  voltage  from  the  anemometer  system. 

TS1  indicates  that  this  correction  is  valid  over  temperature  differences  of  approximately 

50  °C. 


1 .  Turn  on  1FA300,  data  acquisition  PC  and  start  ThermalPro®  software,  select 
Calibration/Probe  Data  from  the  toolbar,  see  Figure  5. 

2.  Connect  appropriate  length  co-axial  cable  and  probe  holder  to  desired  input 
channel  on  the  1FA300  mainframe. 

3.  Install  the  IFA300  thermocouple  in  a  suitable  location  to  measure  the  tunnel 
air  temperature.  This  temperature  will  become  the  calibration  temperature. 

4.  Install  a  shorting  probe  into  probe  holder  and  press  the  "Read"  button  on  the 
Calibration/Probe  Data  window,  see  Figure  A.  1 .  This  will  read  the  cable  and  probe  holder 
resistance. 


5.  Fill  in  the  additional  information  from  the  selected  hot  wire/film  storage  box, 
such  as  sensor  serial  number,  sensor  type,  wire  or  film  and  operating  resistance.  See  the  TS1 
IFA300  operator's  manual  for  a  discussion  on  setting  the  offset  and  gain.  For  fluid  temperature 
reading  select  the  appropriate  channel  that  the  thermocouple  is  plugged  in  to  on  the  back  of  the 
IFA300,  usually  channel  A.  The  calibration  method  is  typically  “Acquire  E  and  Type  Velocity" 
in  a  manual  mode.  Thus  the  1FA300  will  automatically  read  the  output  voltage,  and  the  operator 
must  type  in  the  velocity  for  each  calibration  data  point.  Once  the  probe  data  sheet  is  filled  in. 
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press  the  green  “Calibrate”  button,  which  will  display  the  “Calibration  Conditions  Setup” 
window,  see  Figure  A. 2. 


Figure  AM .  ThermalPro®  software,  Calibration/Probe  Data  Window 


Figure  AI-2.  ThermalPro®  software.  Calibration  Conditions  Setup  Window 
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Figure  AI-3.  Hot  Wire/film  Installed  in  TSI  Calibration  Wind  Tunnel 

6.  Replace  the  shorting  probe  with  the  hot  wire/film  sensor  to  be  calibrated  and 
install  the  probe  holder  and  sensor  into  the  TSI  Calibration  Wind  Tunnel.  Align  the  sensor  in  the 
center  of  the  test  section  in  the  proper  orientation  to  the  flow,  as  shown  in  Figure  A. 3. 

7.  Fill  in  the  additional  information  needed  on  the  “Calibration  Conditions 
Setup”  window  such  as  atmospheric  pressure,  tunnel  air  temperature  and  the  anticipated  min  and 
max  velocity  range  for  which  the  calibration  will  cover,  taking  care  to  select  the  correct  units 

8.  Referencing  the  Hotfilm  Calibration  Spreadsheet  (Table  2  in  main  text) 
configure  the  tunnel  and  pressure  reading  equipment  for  the  first  velocity,  switch  power  onto  the 
tunnel  and  set  the  first  pressure/velocity. 

9.  Manually  input  the  velocity  from  the  spreadsheet  into  the  “Calibration  » 

Conditions  Setup”  window  and  press  the  “Acquire”  button.  The  indicator  next  to  this  button  will 

turn  read  while  the  voltage  is  being  read. 

10.  Repeat  steps  8  and  9  for  each  velocity  in  the  spreadsheet.  When  finished 
press  the  green  “Next  Screen"  button  and  the  acquired  data  will  be  displayed  in  tabulated  form. 

1 1 .  Pressing  the  “Curve”  button  displays  a  graph  of  the  calibration  data,  a  4lh 
order  polynomial  fit  of  the  data,  and  the  %  error  between  the  data  and  the  curve  fit.  An  example 
of  curve  is  shown  in  Figure  A. 4. 
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Fig  AI-4.  Hot  Wire/Film  Calibration  Plot 

Once  calibrated  the  hot  wire  anemometer  is  available  for  use  in  experiments.  A 
drawback  to  the  research  grade  hot  wire/film  sensors  is  that  they  require  extreme  care  in  use. 
They  are  easily  damaged  and/or  broke.  If  this  occurs  then  a  new  sensor  has  to  be  acquired  and 
the  calibration  process  repeated.  For  this  reason,  it  usually  makes  good  sense  to  calibrate  several 
sensors  so  that  they  are  available  in  case  one  is  broken,  and  thus  testing  is  not  interrupted. 
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APPENDIX  II 


VELOCITY  CHARACTERIZATION  PROCEDURE 

A  set  procedure  was  followed  in  measuring  the  velocity  profiles.  The  first  step 
was  to  securely  install  the  wind  tunnel  to  be  characterized  onto  the  laboratory  workbench  to 
reduce  any  motion  between  the  tunnel  and  the  height  gage  (24-in  Mitutoyo  Series  192)  used  to 
position  the  hot  wire  anemometer  above  the  tunnel  floor.  The  height  gage  was  positioned  over 
the  desired  position. 

Determining  the  reference  height  of  the  hot  wire  anemometer  was  an  extremely 
crucial  step  in  the  measurement  process  for  two  reasons.  First,  it  provided  an  accurate  position 
of  the  probe  above  the  wind  tunnel  floor  so  that  height  vs.  velocity  correlations  could  be  made. 
Second,  it  provided  protection  from  inadvertently  damaging  the  probe  by  running  it  into  the 
floor.  To  accomplish  this  zeroing  of  the  height  gage,  an  old  broken  probe  was  modified  with  a 
flat  machined  tip.  The  difference  in  length  from  this  machined  probe  and  the  actual  probe  to  be 
used  for  measurements  was  determined  by  using  an  optical  comparator.  This  device  enlarges 
small  objects  and  permits  measurements  to  within  25  pm.  The  selection  of  common  reference 
points  on  each  probe  body  allowed  the  difference  in  their  lengths  to  be  determined.  The  machine 
probe  was  placed  in  the  probe  holder  on  the  height  gage  and  the  height  gage  was  lowered  until 
the  machined  probe  just  touched  the  piston  on  the  tunnel  floor.  At  this  point  the  height  gage  was 
zeroed  and  the  process  of  raising  and  lowering  the  machined  probe  was  repeated  to  verify  that  a 
reproducible  zero  has  been  obtained.  Once  satisfied  that  the  zero  height  was  reproducible,  the 
height  gage  was  raised  to  the  length  difference  between  the  machined  and  actual  probe 
determined  from  the  optical  comparator.  This  height  was  now  the  true  zero  height  for  the  actual 
probe  and  the  height  gage  is  once  again  zeroed  at  this  position.  Finally  the  height  gage  was 
raised  and  the  machined  probe  removed. 

Air  conditioned  to  the  required  temperature  and  humidity  was  routed  from  the 
Miller-Nelson  ECU  system  to  the  tunnel  inlet.  At  the  same  time,  heating  tape  surrounding  the 
tunnel  was  turned  on  and  the  tunnel  allowed  to  equilibrate  to  the  appropriate  test  conditions. 

This  often  required  several  hours  to  achieve.  When  this  was  accomplished,  the  associated 
measuring  instrumentations  and  computers  are  turned  on  and  allowed  to  warm  up  including  the 
ThermalPro®  system.  From  the  latter’s  menu  bar  “Acquisition”  and  “Probe  Table”  was  selected. 
The  associated  computer  monitor  screen  shown  in  Figure  AIL  1  was  displayed.  After  verifying 
that  the  correct  sensor  is  displayed  in  the  probe  table,  the  cable  and  probe  resistance  were 
measured  by  installing  the  shorting  plug  in  the  probe  holder  and  pressing  the  “Read  Cable” 
button.  Once  read,  the  software  prompted  replacement  of  the  shorting  plug  with  the  appropriate 
sensor.  Then  the  “Read  Probe”  button  was  pressed  and  the  probe  resistance  read.  After 
verifying  that  the  other  settings  listed  are  correct  and  any  values  changed,  the  “SaveLine”  button 
was  pressed  and  the  Probe  Table  was  updated  with  the  current  values.  Finally  the  “Rename” 
button  in  the  experiment  area  of  the  screen  was  pressed  and  a  file  name  entered  for  the  data  files 
saved  during  the  experiment.  When  this  was  completed,  the  “Next  Screen”  button  was  pressed 
to  proceed  to  the  “Acquisition  Conditions  Setup”  screen  as  shown  in  Figure  II. 2. 
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Figure  AIM.  TSI  ThermalPro®  Software;  Acquisition  -  Probe  Table  Screen 


Figure  AII-2.  TSI  ThermalPro®  Software,  Acquisition  Conditions  Setup  Screen 
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The  “Acquisition  Conditions  Setup"  screen  listed  the  settings  that  were  to  be  used 
for  acquiring  the  velocity  profdes.  The  data  on  the  left  side  of  this  screen  was  updated  as 
required.  The  “Acquisition  Control"  area  on  the  right  side  of  the  screen  determined  how  much 
and  how  fast  the  velocity  readings  were  to  be  taken.  Typically,  1,000  readings  over  a  1  second 
time  interval  was  found  to  be  adequate. 

At  this  point  in  the  test  setup,  the  hot  wire  sensor  was  installed  in  the  probe 
holder.  The  sensor  was  lowered  into  the  tunnel  and  the  flexible  rubber  seal  (see  Figure  26) 
placed  between  the  sensor  and  the  test  section  ceiling.  This  seal  prevented  air  from  entering  or 
exiting  the  tunnel  through  the  instrument  port  that  might  otherwise  have  influenced  the  velocity 
profile  measurements.  The  sensor  was  positioned  at  0.5  cm  above  the  wind  tunnel  floor  and  the 
final  flow  adjustments  were  made  prior  to  measuring  the  profile.  The  “Acquisition  Control 
Mode"  was  set  to  “View  Only”  and  the  air  velocity  was  monitored.  The  tunnel  flow  rate  was 
adjusted  until  the  velocity  was  equal  to  the  operational  value  at  the  0.5-cm  height  for  the  selected 
velocity  condition  (low,  medium  or  high).  Once  the  correct  velocity  at  0.5  cm  had  been 
established,  the  various  parameters  such  as  Miller-Nelson  settings,  flow  meter  settings,  tunnel 
humidity,  temperature  and  static  pressure  were  recorded. 

The  “Acquisition  Control  Mode"  was  then  switched  back  to  “Write  Only"  and  the 
hot  wire  sensor  lowered  to  the  lowest  measurement  height  of  1 .27mm  above  the  test  section 
floor.  Although  it  was  desirable  to  go  closer  to  the  surface,  this  limiting  height  was  selected  to 
avoid  a  heat  convection  phenomenon  that  affected  the  hot  wire  output  below  this  height  as  well 
as  to  prevent  inadvertent  damage  to  the  sensors  if  it  came  in  contact  with  the  tunnel  floor.  A 
profile  reading  was  obtained  at  each  of  the  2 1  probe  heights  covering  a  range  of  0. 1 27  to  4.5  cm. 
The  complete  listing  of  probe  heights  is  shown  in  Table  II.  1 . 

Table  All- 1 .  Probe  Heights  for  5-cm  Wind  Tunnel  Characterization 


Probe  Height  (y) 

(in) 

(cm) 

0.050 

0.127 

0.098 

0.25 

0.148 

0.376 

0.197 

0.5 

0.246 

0.625 

0.295 

0.75 

0.394 

1 

0.492 

1.25 

0.591 

1.5 

0.689 

1.75 

Probe  Height  (y) 

0.787 

2 

0.886 

2.25 

0.984 

2.5 

1.083 

2.75 

1.181 

3 

1.280 

3.25 

1.378 

3.5 

1.476 

3.75 

1.575 

4 

1.673 

4.25 

1.772 

4.5 

After  the  IFA300  recorded  a  measurement,  the  probe  was  raised  to  the  next 
height.  This  was  repeated  until  the  last  height  of  4.5  cm.  At  this  point  the  measurement  process 
was  over  and  the  stored  readings  were  analyzed  by  the  Thermal  Pro®  software.  This  involved  a 
number  of  statistical  calculations.  However,  the  two  most  critical  results  were  the  average  total 
velocity  and  the  turbulence  intensity.  The  final  step  was  to  build  the  velocity  field,  in  which  all 
of  the  readings,  which  up  to  this  point  have  been  stored  in  a  number  of  separate  files,  were 
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combined  into  tabulated  and  plotted  formats.  Typically,  the  Thermal  Pro®  plots  were  not  used 
and  the  tabulated  data  was  imported  into  Microsoft  Excel®  for  comparison  plotting  with  other 
profiles.  Normally,  two  plots  were  made  at  each  position  and  wind  speed  setting.  The  hot  wire 
anemometer  was  then  repositioned  to  the  0.5-cm  height  and  the  process  repeated. 

Care  must  be  taken  in  interpreting  the  hot  wire  velocity  measurements  as  the 
approaches  the  floor.  Heat  transfer  effects  due  to  the  proximity  of  the  metal  floor  wall  can 
the  temperature  of  the  hot  wire  giving  false  velocity  readings.  The  closest  measurement 
test  section  floor  for  the  5-cm  Wind  Tunnel  characterization  was  0.050  in.  ( 1 .27  mm). 


probe 
affect 
to  the 
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APPENDIX  III 


TABULATED  VELOCITY  PROFILES 


LIST  OF  FIGURES 

Figure  AIII-1 .  5-cm  Tunnel  Operational  Profiles 
Figure  All  1-2.  Tunnel  3A  (Low  and  Medium  Velocities) 

Figure  AIII-3.  Tunnel  3B  (Low  Velocity) 

Figure  AIII-3.  (Continued).  Tunnel  3B  (Low  Velocity) 

Figure  AIII-3.  (Continued).  Tunnel  3B  (Medium  Velocity) 

Figure  AIII-3.  Continued).  Tunnel  3B  (Medium  Velocity  Continued) 

Figure  AIII-3.  (Continued).  Tunnel  3B  (High  Velocity) 

Figure  AIII-3.  (Continued).  Tunnel  3B  (High  Velocity  Continued) 

Figure  AIII-4.  Tunnel  3C 
Figure  AIII-5.  Tunnel  3D 
Figure  AIII-6.  Tunnel  3E 
Figure  AIII-7.  Tunnel  3F  (Low  Velocity) 

Figure  AIII-7.  (Continued).  Tunnel  3F  (Low  Velocity  Continued) 

Figure  AIII.7.  (Continued).  Tunnel  3F  (Medium  Velocity) 

Figure  AIII.7.  (Continued).  Tunnel  3F  (Medium  Velocity  Continued) 

Figure  AIII.7.  (Continued).  Tunnel  3F  (High  Velocity) 

Figure  All  1.8.  Tunnel  3F  -  Effect  of  Temperature 

Figure  All  1.9.  Tunnel  3F  -  Effect  of  Surface  Roughness  (Low  Velocity) 

Figure  AIII.9.  (Continued).  Tunnel  3F  -  Effect  of  Surface  Roughness  (Medium  Velocity) 

Figure  AIII.9.  (Continued).  Tunnel  3F  -  Effect  of  Surface  Roughness  (High  Velocity) 

Figure  AIII-10.  Tunnel  3G 

Figure  AIII-10.  (Continued).  Tunnel  3G 

Figure  AIII-1 1.  Tunnel  31 

Figure  AIII-12.  Tunnel  3J 

Figure  AIII-1 3.  Tunnel  3K 

Figure  AIII-14.  Tunnel  3L  (Low  Velocity) 

Figure  AIII-14.  (Continued).  Tunnel  3L  (Medium  Velocity) 

Figure  AIII-14.  (Continued).  Tunnel  3L  (High  Velocity) 
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Operational  Velocity 

Probe 

Height 

Low 

Medium 

High 

(cm) 

m/s 

m/s 

m/s 

0.010 

0.003 

0.074 

0.265 

0.016 

0.004 

0.121 

0.434 

0.027 

0.007 

0.198 

0.713 

0.044 

0.012 

0.326 

1.146 

0.072 

0.020 

0.534 

1.609 

0.119 

0.033 

0.778 

2.071 

0.195 

0.053 

1.022 

2.533 

0.320 

0.173 

1.325 

2.825 

0.500 

0.194 

1.431 

3.026 

0.525 

0.196 

1.446 

3.054 

0.862 

0.218 

1.557 

3.264 

1.000 

0.222 

1.580 

3.308 

1.414 

0.236 

1.649 

3.439 

2.321 

0.253 

1.740 

3.611 

Figure  A 1 1 1- 1 . 
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Tunnel  3A 

Position  A 

Q130T35H50R2 

Turb  Int. 

(%) 

17.964 

16.425 

Is- 

CO 

CO 

269  31- 

11.523 

9.705 

8.023 

6.934 

5.684 

5.418 

4.739 

4.485 

3.505 

3.556 

3.612 

4.069 

3.901 

CO 

N- 

o 

in 

8.012 

10.435 

16.770 

Velocity 

(m/s) 

9060 

1.174 

1.295 

1.364 

CO 

CO 

1.469 

1.546 

Is- 

CO 

CO 

1.651 

1.713 

1.726 

1.742 

1.740 

1.729 

1.738 

1.720 

069'  t 

1.672 

00 

oo 

in 

1.444 

1.140 

Q130T35H50R1 

Turb  Int. 

(%) 

19.044 

14.874 

15.050 

13.059 

10.688 

10.514 

8.947 

7.322 

5.773 

5.307 

5.385 

CO 

CO 

CO 

4.052 

3.954 

O 

CO 

CO 

CO 

3.312 

3.602 

5.233 

7.496 

10.778 

15.033 

Velocity 

(m/s) 

8060 

CO 

CNJ 

co 

o 

CO 

CO 

1.471 

1.526 

1.582 

1.679 

1.738 

1.732 

1.759 

CO 

oo 

h- 

CM 

00 

1^ 

1.765 

1.765 

1.735 

1.756 

1.685 

CD 

5 

CO 

r^ 

1.179 

Tunnel  3A 

Position  A 

Q11T35H50R2 

Turb  Int. 

2.936 

4.557 

5.693 

5.215 

5.667 

Z88S 

6.142 

5.736 

00 

CO 

5.624 

6.213 

5.248 

5.976 

688  9 

5.097 

o 

CO 

00 

4.737 

4.177 

4.137 

3.577 

3.622 

Velocity 

(m/s) 

CO 

CO 

o 

0.145 

in 

CO 

d 

CO 

00 

d 

o 

o 

CM 

d 

0.212 

0.220 

o 

CNJ 

CNJ 

d 

0.209 

0.197 

0.186 

0.174 

co 

CO 

d 

0.162 

o 

CO 

d 

0.154 

0.150 

0.147 

0.143 

0.137 

IZl  0 

Q11T35H50R1 

Turb  Int. 

(%) 

2.831 

4.025 

CO 

4.258 

4.704 

CO 

220  9 

5.249 

CO 

1^ 

4.527 

4.447 

5.430 

co 

co 

o 

in 

5.082 

5.005 

5.640 

4.992 

4.719 

4.822 

4.174 

3.029 

Velocity 

(m/s) 

0.131 

CO 

d 

0.169 

CO 

co 

d 

0.203 

0.213 

0.222 

022  0 

0.210 

0.200 

0.187 

0.177 

0.169 

0.162 

0.157 

0.153 

0.150 

0.147 

CO 

d 

1^ 

CO 

d 

0.123 

C/i 

<L> 


CJ 

JO 

<L> 

> 

E 

2 

-5 
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TD 

C 

03 
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< 
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13 

c 

a 

3 

H 

(N 


< 

0> 
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D 
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in 
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in 
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CO 
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-4— » 

sz 

U) 

E 

Q 

CNJ 
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i^ 

co 

d 

CM 
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h- 

CM 
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CM 

co 

co 

h- 

CO 
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'cd 
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o 

co 

co 

i^ 

co 

in 

CM 

CD 
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^r 

co 

o 

co 

co 

m 

co 

CM 
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in 

CD 

CD 

CD 

CD 

CD 

CD 

co 

oo 

oo 

CO 

00 

co 

oo 

Is- 
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Is- 

1^ 

2 
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o 

o 

T— 

CM 

CM 

CO 

in 
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r- 
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CM 

CL 

O 

LO 

m 

Turb  Int. 

1.703 

2.696 

3.286 

4.035 

4.687 

5.078 

CO 

GO 

^r 

4.901 

4.403 

4.291 

4.981 

4.242 

OO 

CO 

4.690 

4.213 

3.702 

4.170 

3.411 

3.636 

3.239 

2.847 

QQ 

CO 

~Q> 

CD 

a 

o 

CO 

1 — 

CM 

O 

Velocity 

(m/s) 

0.149 

0.149 

0.160 

0.177 

0.190 

0.204 

0.230 

0.239 

0.238 

0.228 

0.216 

0.204 

0.194 

0.186 

0.183 

0.178 

0.176 

0.174 

0.171 

0.166 

0.156 

C 

c 

1— 

c/d 

o 

CL 

CL 

CD 

LO 

m 

LO 

Turb  Int. 

1.268 

2.220 

3.044 

oo 

co 

CO 

3.955 

3.955 

4.034 

4.014 

4.040 

3.592 

3.934 

4.497 

4.417 

4.099 

4.113 

4.259 

CO 

4.097 

3.998 

3.688 

2.589 

h— 

CM 

O 

Velocity 

(m/s) 

0.150 

0.148 

0.160 

0.173 

0.191 

0.205 

0.225 

0.238 

0.240 

0.228 

0.218 

0.206 

0.192 

0.186 

0.181 

0.177 

0.175 

0.174 

0.171 

0.166 

0.155 

CM 

CL 

o 

LO 

hi 

LO 

Turb  Int. 

o'" 

2.522 

3.316 

CD 

CO 

4.616 

4.883 

5.478 

4.950 

4.929 

5.439 

4.698 

5.631 

4.637 

4.725 

4.511 

4.069 

4.525 

4.300 

4.150 

3.666 

3.855 

OO 

OO 

r— 

CM 

CD 

CO 

“a> 

c 

o 

C 

O 

1 - 

CM 

O 

Velocity 

(m/s) 

0.131 

0.136 

0.153 

0.170 

0.188 

0.201 

0.222 

0.228 

0.221 

0.216 

0.207 

0.196 

0.183 

0.174 

0.167 

0.161 

0.158 

0.152 

0.149 

0.143 

0.136 

c 

13 

I— 

co 

o 

CL 

cl 

CO 

LO 

n= 

Turb  Int. 

2.195 

3.195 

4.287 

4.501 

5.106 

4.351 

Is- 

oo 

CO 

4.194 

4.145 

4.614 

4.802 

4.375 

OO 

CO 

or 

or 

4.462 

4.716 

4.656 

4.631 

CO 

CM 

OO 

•^r 

r^- 

CD 

OO 

CO 

4.042 

2.930 

LO 

CO 

I— 

CM 

O 

Velocity 

(m/s) 

0.129 

0.137 

0.154 

0.170 

0.187 

0.205 

0.218 

0.225 

0.224 

0.217 

0.206 

OO 

CD 

CD 

0.186 

0.175 

0.171 

0.164 

0.160 

0.155 

0.150 

0.146 

0.137 

CO 

CL 

CO 

LO 

HI 

Turb  Int. 

CT" 

2.945 

4.246 

4.408 

5.496 

5.420 

5.795 

5.683 

5.174 

5.705 

5.729 

6.121 

5.472 

5.633 

5.255 

4.515 

4.474 

4.161 

3.947 

4.007 

3.682 

3.209 

CD 

CO 

"a> 

< 

c 

o 

CO 

1 — 
CM 

O 

Velocity 

(m/s) 

0.126 

0.139 

0.161 

0.179 

0.197 

0.207 

0.220 

0.218 

0.205 

0.194 

0.181 

0.170 

0.161 

0.153 

0.147 

0.143 

0.140 

0.137 

0.139 

0.139 

0.133 

£= 

C 

3 

K— 

<75 

O 

CL 

CM 

CL 

CO 

LO 

HI 

Turb  Int. 

o'" 

2.152 

3.261 

CM 

CO 

CO 

CO 

OO 

CD 

CO 

4.340 

4.144 

4.280 

4.297 

3.888 

4.202 

3.971 

4.053 

4.314 

3.975 

3.401 

3.327 

3.046 

3.478 

2.944 

2.714 

2.495 

CO 

1— 

CM 

O 

Velocity  1 

(m/s) 

0.128 

0.142 

0.163 

CO 

OO 

CO 

0.197 

0.208 

0.222 

0.217 

0.210 

0.194 

0.179 

0.166 

0.156 

0.147 

0.140 

0.136 

0.132 

0.130 

0.131 1 

0.130 

0.124 

(cm) 

0.127 

LO 

CM 

CD 

0.376 

LO 

CD 

0.625 

LO 

Is- 

CD 

LO 

CM 

LO 

LO 

Is— 

CM 

LO 

CM 

CM 

LO 

CM 

LO 

r— 

CM 

CO 

LO 

CM 

CO 

LO 

OO 

LO 
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CO 

■^r 

LO 

CM 
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CD 

T — 
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CO 
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LO 
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CD 

CD 
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CO 
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LO 
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CD 

CD 
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CD 
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Tunnel  3B 

Position  E 

Turb  Int. 

_ (%) _ 

o 

CD 

CD 

CN 

4.348| 

5.123 

5.519 

6989 

5.997 

7.025 

6.433 

6.929 

5.930 

6.324 

6.124 

4.850 

4.766 

4.597 

CN 

O' 

O; 

O' 

4.797 

5.014 

4.919 

4.281 

3.730 

Velocity 

CO 

E, 

CD 

O' 

O 

0.154 

o 

CD 

00 

d 

0.205 

0.218 

0.235 

0.242 

Is- 

co 

CN 

d 

0.225 

0.216 

0.204 

0.197 

0.189 

0.182 

CD 

r- 

d 

0.172 

0.170 

C30 

CD 

d 

0.167 

0.158 

Turb  Int. 

g 

05 

CNJ 

CD 

CNJ 

3 

00 

O' 

i^ 

O' 

o 

ID 

CD 

LD 

Is- 

Is- 

LD 

LD 

O' 

05 

LD 

T— 

ID 

ID 

CO 

O' 

LD 

00 

O' 

oq 

O' 

O' 

O' 

CO 

ID 

Is- 

o 

LD 

ID 

Is- 

CO 

LD 

ID 

00 

05 

O 

LD 

O 

00 

05 

ID 

Is- 

o 

CD 

ID 

05 

O' 

LD 

CD 

CN 

O' 

LD 

ID 

O 

LD 

ID 

CO 

r- 

Is- 

ID 

CN 

O' 

00 

O' 

o 

05 

05 

CO 

Velocity 

"c/T 

E, 

ID 

CO 

O 

O' 

O' 

o 

ID 

CD 

T— 

d 

S 

T — 

o 

05 

05 

T— 

d 

LD 

CNJ 

d 

00 

CNJ 

CNJ 

d 

O' 

CNJ 

d 

Is- 

co 

CNJ 

d 

0.225 

O' 

CN 

d 

0.202 

1 

d 

0.187 

5 

d 

0.176 

CN 

Is- 

d 

0.168 

0.167 

O' 

CD 

d 

0.156 

Tunnel  3B 

Position  D 

Q12T35H50R2 

Turb  Int. 

Co' 

ID 

ID 

O 

CO 

o 

CD 

CD 

O' 

ID 

CO 

O 

CD 

CO 

ID 

05 

LD 

o* 

CD 

CO 

CD 

00 

CO 

O' 

d 

00 

00 

CD 

CNJ 

Is- 

CO 

CD 

ID 

CO 

CD 

00 

00 

05 

CD 

o 

LD 

O 

CD 

o 

05 

00 

LD 

CN 

00 

05 

LD 

oo 

r- 

CD 

LD 

O' 

T— 

LD 

4.798 

4.599 

4.175 

s 

CO 

O' 

o 

CD 

o1 

o 

CD 

Is- 

co 

Velocity 

'in 

i 

LO 

O 

r- 

CD 

o 

O 

05 

O 

Is- 

O 

CNJ 

d 

o- 

CNJ 

CNJ 

d 

CD 

CO 

CNJ 

d 

LD 

O' 

CNJ 

d 

00 

CO 

CN 

d 

CN 

CN 

CN 

d 

O' 

O 

CN 

d 

00 

00 

d 

O' 

1^ 

d 

CO 

CD 

d 

ID 

LD 

d 

05 

O' 

d 

O' 

O' 

d 

CO 

O' 

T— 

d 

CO 

O' 

T— 

d 

O' 

O' 

T— 

d 

CD 

O' 

d 

CO 

O' 

d 

Q12T35H50R1 

Turb  Int. 

g 

CO 

CO 

4.506 

5.516 

o- 

CNJ 

CD 

LD 

6.214 

CD 

o 

o 

CD 

5.302 

CNJ 

05 

ID 

CN 

LD 

4.573 

s 

o 

CD 

05 

CD 

O' 

LD 

O' 

05 

05 

ID 

5.780 

CO 

CO 

CD 

LD 

O' 

CD 

05 

O' 

ID 

r- 

LD 

ID 

C30 

CN 

ID 

Is- 

i^ 

CN 

LD 

4.516 

CD 

ID 

CO 

Velocity 

in 

E_ 

o 

ID 

O 

Is- 

CD 

T— 

o 

05 

T— 

o 

h- 

O 

CNJ 

d 

1^ 

CNJ 

CNJ 

d 

00 

CO 

CNJ 

d 

CO 

O' 

CNJ 

d 

h- 

co 

CNJ 

d 

CNJ 

CN 

d 

CO 

o 

CN 

d 

00 

00 

d 

O' 

1^ 

d 

O' 

CD 

d 

ID 

ID 

d 

05 

O' 

d 

ID 

O' 

d 

CO 
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O' 

d 

Is- 

o- 

d 

CO 

O' 

d 

1^ 

ID 

CD 

ID 

ID 

ID 

ID 

LD 

LD 

CN 

LD 

LD 

LD 

CO 

ID 

ID 

ID 

O' 

ID 

ID 

^  CN 

E  ^ 

CN 

d 

Is- 

co 

d 

CN 

CD 

1^ 

d 

CN 

1^- 

CN 

CN 

CN 

h- 

CN 

CN 

CO 

CO 

Is- 

CO 

CN 

O' 

O' 

o 

d 

d 

o 

00 

00 

1^ 

CD 

LD 

O' 

CN 

05 

1^ 

CD 

O' 

CO 

O 

00 

CD 

ID 

CO 

CN 

ID 

05 

O' 

05 

O' 

05 

05 

05 

05 

00 

00 

00 

00 

00 

00 

00 

1^ 

Is- 

r- 

N> 

r— 

o 

T— 

CN 

CN 

CO 

O' 

LD 

CD 

1^ 

00 

05 

O 

T— 

CN 

CO 

O; 

LD 

p 

Is- 

£  d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

T— 

T 

T 

T— 

APPENDIX  III 


93 


Tunnel  3B 

Position  B 

CM 

o' 

o 

m 

X 

in 

Turb  Int. 

(%) 

17.497 

16.829 

15.784 

13.479 

10.335 

9.394 

8.275 

8.397 

9.007 

9.424 

10.517 

9.637 

CD 

11.689 

9.661 

CD 

OO 

OO 

r- 

6908 

6.497 

9.138 

12.580 

19.912 

CO 

1 — 
CO 

O 

Velocity 

(m/s) 

609  0 

O 

CD 

O 

1.252 

1.336 

1.394 

1.450 

1.528 

1.508 

1.490 

1.512 

1.449 

CD 

1.459 

1.505 

1.537 

1.599 

O 

CD 

1.598 

1.510 

1.399 

1.084 

Q143T35H50R1 

Turb  Int. 

18.113 

16.431 

14.419 

11.383 

10.590 

9.815 

8.952 

7.779 

8.927 

11.129 

oo 

o 

o 

— 

9.927 

O 

o 

OO 

o 

9.353 

8.711 

6062 

7.310 

9689 

8.043 

12.156 

19.241 

Velocity 

(m/s) 

0.822 

1.123 

1.266 

1.342 

1.424 

1.439 

O 

OO 

1.512 

1.536 

1.471 

1.509 

1.484 

1.491 

1.520 

1.598 

1.585 

1.594 

1.594 

1.544 

O 

co 

1.129 

Tunnel  3B 

Position  C 

Q143T35H50R2 

Turb  Int. 

(%) 

18.866 

14.917 

13.079 

11.817 

10.824 

9.620 

6998 

7.373 

7.526 

3 

oo 

6.855 

OO 

o 

oo 

CD 

6.763 

6.216 

1^- 

CD 

O 

CD 

6.205 

6.725 

6.851 

7.597 

9.673 

16.323 

Velocity 

(m/s) 

0.907 

CD 

1.315 

1.361 

CD 

CO 

m 

oo 

1.547 

1.566 

1.611 

Z60L 

CM 

CO 

CD 

CD 

CD 

1.647 

1.655 

1.646 

1.643 

1.646 

O 

O 

CD 

1.537 

1.420 

1.131 

Q143T35H50R1 

Turb  Int. 

(%) 

17.594 

15.543 

13.042 

10.609 

10.655 

8.538 

9.144 

699  L 

OO 

CD 

7.216 

8.032 

7.322 

7.134 

7.513 

OO 

o 

CD 

6.240 

6.811 

7.021 

LLVL 

11.829 

17.303 

Velocity 

(m/s) 

0.962 

1.193 

1.322 

CD 

O 

1.419 

1.509 

1.546 

1.616 

1.621 

1.662 

1.618 

1.661 

OO 

CD 

CD 

1.644 

1.649 

1.619 

1.627 

1.625 

1.518 

1.410 

1.169 

Tunnel  3B 

Position  A 

Q143T35H50R3 

Turb  Int. 

(%) 

18.277 

14.201 

12.412 

10.422 

9.251 

8.149 

7.219 

6.442 

5.960 

4.989 

4.395 

4.157 

4.061 

3.224 

CD 

OO 

CO 

3.611 

3.873 

4.845 

7.041 

10.200 

16.420 

Velocity 

(m/s) 

0.975 

1.267 

1.342 

1.458 

1.496 

1.535 

CO 

oo 

in 

1.650 

1.666 

1.701 

1.702 

oo 

1.714 

1.738 

OO 

o 

1.702 

1.692 

1.643 

1.605 

1.495 

1.201 

Q143T35H50R2 

Turb  Int. 

c %: r 

17.242 

14.154 

CD 

O 

CO 

10.915 

10.055 

8.519 

7.378 

6.242 

5.329 

5.035 

4.398 

4.558 

3.979 

3.812 

3.344 

3.774 

4.119 

5.168 

6.657 

10.654 

16.871 

Velocity 

(m/s) 

1.029 

1.266 

CM 

o 

1.441 

1.495 

1.562 

1.606 

1.640 

1.668 

1.724 

1.711 

1.724 

1.714 

1.724 

1.724 

1.709 

1.687 

1.685 

1.613 

OO 

1.227 

m 

CD 

m 

m 

m 

in 

in 

in 

CM 

m 

in 

m 

CO 

in 

in 

m 

in 
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h- 

o 
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r— 
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T — 
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CM 
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00 
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r-~ 
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CD 
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oo 
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in 

CO 

CM 

p 

CD 

CD 

OT 

CD 

CD 

CD 
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oo 
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Figure  AIII-3.  (Continued).  Tunnel  3B  (Medium  Velocity  Continued) 
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Figure  AIII-3.  (Continued).  Tunnel  3B  (High  Velocity  Continued) 
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Tunnel  3F 

Position  E 

Q16T35H40R2 

Turb  Int. 

(%) 

o 

o 

CO 

o 

1.053 

oo 

o 

CO 

CM 

1.434 

1.463 

1.492 

1.616 

1.621 

CO 

LO 

o 

CO 

1.521 

1.370 

o 

Is- 

CM 

1.256 

1.217 

1.170 

1.127 

oo 

o 

Z960 

r- 

o 

d 

Velocity 

c/> 

i 

0.145 

0.162 

CO 

00 

d 

O 

CM 

d 

CM 

d 

0.223 

0.228 

022  0 

Is- 

o 

CM 

d 

0.194 

CO 

00 

d 

CM 

r- 

d 

0.165 

O 

CO 

d 

0.157 

0.155 

CO 

in 

7 — 

o 

CM 

in 

d 

in 

d 

00 

d 

d 

Q16T35H40R1 

Turb  Int. 

vO 

0.572 

CO 

CM 

CD 

d 

co 

CO 

CM 

CM 

CD 

00 

CM 

o 

CD 

CM 

Is- 

CO 

CO 

CD 

CO 

^r 

CO 

co 

s 

in 

in 

co 

CO 

00 

CO 

CO 

CO 

1^ 

CM 

CO 

CM 

00 

00 

CD 

co 

co 

00 

o 

0.949 

CD 

CO 

d 

Velocity 

</) 

E, 

CO 

XT 

7 — 

o 

co 

d 

00 

d 

2020 

CO 

CM 

d 

CM 

CM 

d 

0.229 

CM 

CM 

d 

0.209 

0.195 

CM 

00 

d 

0.172 

0.165 

0.160 

Is- 

in 

d 

0.155 

0.154 

co 

in 

d 

CM 

in 

d 

0.149 

CM 

d 

Tunnel  3F  _ 

Position  D 

Q16T35H40R2 

Turb  Int. 

2 

o 

CD 

CO 

d 

CD 

d 

r- 

o 

CO 

CD 

o 

CM 

00 

Is- 

CO 

CM 

o 

CO 

00 

^r 

CO 

in 

CD 

CO 

r- 

CD 

O 

CO 

CO 

CM 

CM 

CM 

CM 

8660 

es60 

0.907 

0.748 

0.545 

Velocity 

if) 

1 

0.210 

0.240 

0.272 

0.296 

CO 

d 

0.319 

0.319 

CO 

o 

CO 

d 

CM 

CD 

CM 

d 

0.273 

in 

in 

CM 

d 

00 

co 

CM 

d 

CM 

CM 

d 

CM 

d 

CO 

O 

CM 
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o 

CM 
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00 

CD 

d 

CO 

CD 

d 

^r 

CD 

d 

CM 

CD 

d 

00 

d 

Q16T35H40R1 

Turb  Int. 

2 

CO 

co 

d 

CO 

co 

°0 
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CD 

d 

CM 
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CD 
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r- 

CM 

o 

CD 

O 

CO 

oo 

CM 

CD 

00 

CO 

CD 

CO 

00 

in 
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CM 
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oo 
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CM 

CO 

086  0 

^r 

CO 

CD 
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CD 

oo 
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CO 

00 
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O 

Is- 

d 

CM 

CO 

in 

d 

Velocity 

in 

1 
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0.242 

0.274 

Is- 

CD 

CM 
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co 
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CM 

CM 

CO 
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CM 

CO 
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CD 
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CO 
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CO 

CD 

CM 
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in  in 

in  cm 

in 

in 

in 

co  in 

in 

in  ^  in 

in 

'P  ™ 

CM 

Is- 

d 

CM 

Is- 

CM  T_ 

r- 

CM 

CM 

CM 

co 

^  N 

t  7— 

d 

00 

CO 

d 

X— 

7 — 

CM 

CM 

00 

co 

H  d 

d 

d 

X 

<1) 

-Q 

O 

i— 

CL 


o 

00 

oo 

r- 

co 

in 

XT 

CM 

■*— 

CD 

Is- 

co 

co 

T— 

o 

00 

CO 

in 

CO 

CM 

in 

CD 

CD 

CD 

CD 

CD 

CD 

oo 

oo 

oo 

00 

00 

00 

00 

Is- 

Is- 

h- 

Is- 

c 

o 

o 

7 — 

CM 

CM 

CO 

in 

CO 

Is- 

00 

CD 

o 

7— 

CM 

CO 

in 

co 

Is- 

d 

d 

o 

d 

d 

d 

d 

d 

d 

d 

d 

d 

d 

7 — 

7- 

7 

7- 

7- 

7- 

7 — 

7- 

-o 

<u 

zs 

c 

C 

o 

U 

G 

o 


<L> 

> 

£ 


m 

C 

5 

3 


-o 

o 


c 

o 

U 


< 

<D 

u- 

3 

5XJ 

Lu 


% 


APPENDIX  III 


102 


i 


Tunnel  3F 

Position  C 

Q137T35H50R2 

Turb  Int. 

(%) 

16.935 

14.284 

CO 

CO 

o 

1 1 .263 

9996 

CD 

CO 

CO 

7.872 

7.379 

7.007 

7.193 

6.541 

o 

N- 

O 

1^ 

6.516 

5.574 

5.927 

6.410 

6.501 

6.511 

7.376 

10.526 

CO 

00 

o 

Velocity 

(m/s) 

966  0 

1.288 

1.395 

1.473 

1.538 

1.573 

1.631 

1.661 

o 

o 

N- 

1.703 

1.733 

1.726 

1.731 

1.756 

1.741 

1.731 

1.727 

1.697 

1.621 

1.525 

1.294 

Q137T35H50R1 

Turb  Int. 

(%) 

o 

CO 

CD 

CD 

15.534 

12.214 

11.615 

0066 

8.714 

co 

7.237 

00 

GO 

CD 

7.506 

7.067 

CM 

CO 

CD 

5.846 

5.601 

6.064 

9869 

6.580 

6.529 

7.390 

10.294 

13.433 

Velocity 

(m/s) 

O 

O 

o 

CO 

CM 

1.448 

00 

CO 

1.563 

989  U 

1.632 

CO 

CD 

1.707 

CO 

1.737 

1.728 

1.790 

1.790 

1.742 

co 

CD 

1.725 

o 

N- 

CD 

1.525 

1.284 

Tunnel  3F 

Position  B 

Q153T35H50R2 

Turb  Int. 

(%) 

15.341 

14.805 

12.560 

9.150 

8.939 

7.779 

7.592 

8.372 

9.162 

9.194 

9.652 

9.419 

9.463 

8.596 

7.871 

6.190 

co 

CD 

7.320 

9.543 

12.764 

20.399 

Velocity 

(m/s) 

0680 

1.137 

1.264 

1.368 

1.410 

1.455 

1.457 

1.456 

1.424 

1.397 

1.410 

1.417 

1.465 

co 

CO 

LO 

1.553 

1.528 

1.506 

1.431 

1.279 

686  0 

Q153T35H50R1 

Turb  Int. 

(%) 

in 

o 

o 

CO 

12.191 

6966 

8.439 

CO 

CO 

CD 

1^ 

00 

co 

CD 

1^ 

00 

1^ 

8996 

8.816 

9996 

9.436 

9.321 

9.155 

8969 

6.774 

co 

CD 

CD 

8699 

9998 

13.470 

o 

d 

CM 

Velocity 

(m/s) 

1.122 

1.276 

1.367 
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Figure  AIM. 7.  (Continued).  Tunnel  3F  (High  Velocity  Continued) 
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Figure  All  1.8.  Tunnel  3F  Measurements  -  Effect  of  Temperature 
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Figure  AIII.9.  Tunnel  3F  -  Effect  of  Surface  Roughness  (Low  Velocity) 
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Tunnel  3K 

Position  A 

Q407T35H23R2 

Turb  Int. 

(%) 

16.077 

12.198 

9.984 

9.193 

8.468 

CO 

CO 

00 

7.210 

5.853 

5.112 

00 

00 

00 

3.091 

2.445 

2.238 

2.421 

CD 

CD 

CD 

csi 

3.257 

3.990 

4.956 

00 

5 

CD 

9.308 

14.813 

Velocity 

(m/s) 

2.399 

2.736 

2.939 

3.039 

3.125 

3.203 

3.324 

3.481 

3.577 

3.652 

3.682 

9698 

3.699 

00 

CD 

CO 

3.649 

3.604 

3.540 

3.434 

3.288 

1^- 

00 

o 

00 

CO 

CO 

CN 

Q407T35H23R1 

Turb  Int. 

(%) 

15.816 

12.267 

CO 

CN 

O 

O 

9.075 

8.704 

8.026 

7.147 

6.209 

4.962 

CO 

o 

2.823 

2.386 

2.288 

2.356 

2.824 

3.477 

4.150 

5.002 

6.833 

9.651 

14.619 

Velocity 

(m/s) 

00fr'2 

5 

h- 

Csi 

2.947 

o 

CO 

o 

CO 

3.112 

3.204 

3.348 

o 

CO 

CO 

3.588 

3.653 

3.697 

3.711 

3.707 

3.691 

9998 

3.598 

3.534 

CD 

CO 

CO 

3.293 

3.112 

2.762 

Tunnel  3K 

Position  A 

Q177T35H32R2 

Turb  Int. 

m 

17.279 

15.214 

13.267 

9.738 

o 

CO 

oo 

d 

0908 

7.628 

5.939 

4.642 

3.620 

2.728 

2.362 

2.202 

ZZZ'Z 

2.561 

3.109 

3.844 

5.100 

8.002 

CO 

CO 

00 

d 

16.318 

Velocity 

(m/s) 

8Z60 

1.256 

1.373 

1.449 

1.469 

1.558 

1.612 

o 

CO 

CO 

1.703 

1.730 

CD 

1^ 

1.745 

1.743 

1.732 

1.719 

o 

o 

1^ 

o 

h- 

CD 

1.627 

1.536 

CD 

CO 

1.024 

Q177T35H32R1 

Turb  Int. 

(%) 

19.445 

16.562 

12.444 

10.858 

CN 

CO 

CO 

d 

8.748 

7.481 

5.564 

4.791 

3.824 

3.011 

2.356 

2.243 

2.142 

ZLVZ 

2.917 

3.516 

5.916 

7.404 

10.905 

CD 

CN 

00 

Velocity 

(m/s) 

0.964 

1.215 

1.358 

CO 

00 

1.542 

1.596 

00 

CO 

CO 

1.698 

1.719 

1.739 

1.744 

1.740 

1.731 

1.722 

1.697 

1.672 

1.612 

00 

s 

CO 

CD 

CO 

CD 

O 

O 

Tunnel  3K 

Position  A 

Q11T35H8R2 

Turb  Int. 

(%) 

9890 

969  0 

o 

o 

CO 

d 

00 

CO 

CN 

d 

0.435 

889 ‘0 

CO 

d 

CO 

d 

0.415 

0.644 

0.502 

0.679 

00 

00 

CO 

d 

CD 

d 

9890 

0.333 

0.543 

0.528 

0.714 

0.297 

oo 

oo 

d 

Velocity 

(m/s) 

0.131 

0.151 

0.175 

0.194 

0.204 

0.210 

0.212 

0.206 

0.200 

0.194 

h- 

00 

d 

0.178 

0.171 

CD 

o 

0.155 

00 

d 

0.141 

0.135 

0.128 

0.120 

o 

d 

Q11T35H8R1 

Turb  Int. 

(%) 

0.641 

CN 

N" 

CO 

O 

0.499 

0.517 

00 

CO 

CN 

d 

0.429 

0.543 

0.428 

0.427 

0Z90 

0.453 

CD 

CD 

d 

8090 

0.246 

0.744 

0.734 

0.541 

0.812 

0.832 

0.442 

0.514 

Velocity 

(m/s) 

0.133 

0.152 

0.175 

0.193 

0.204 

0.210 

0.212 

0.207 

0.200 

0.194 

0.187 

0.179 
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0.165 

0.156 
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Tunnel  3L 

Position  E 

CL 

o 

CM 

X 

c 

-Q 

D 

1- 

I — m 

12.0601 

9.948| 

8.562| 

8.057| 

7.820| 

o 

CD 

r- 

5.937 1 

CO 

CD 

3.683| 

2.817| 

2.202| 

2.058| 

2.031 1 

2.096| 

2.21 4 1 

2.441| 

CD 

CM 

CO 

4.682| 

i^- 

CM 

10.6741 

16.931 1 

CO 

h— 

CO 

CD 

CO 

O 

Velocity 

(m/s) 

2.61 4 1 

2.842| 

in 

CD 

CD 

CM 

3.058] 

o 

o 

CO 

3.201| 

3.336] 

3.433| 

3.472| 

3.482| 

3.486] 

CD 

00 

CO 

3.492| 

3.496] 

3.502] 

3.499] 

3.457] 

CO 

3.295] 

3.082] 

CD 

CO 

CD 

CM 

Position  D 

Q398T35H20R1 

Turb  Int. 

i  i 

14.971 1 

12.268| 

10.239] 

8.733] 

8.163| 

7.504| 

6.281  ] 

4.7601 

o 

CD 

CO 

2.870] 

2.268] 

\zooz 

o 

CM 

O 

CM 

1.924] 

2.083] 

2.190] 

2.972] 

4.084] 

6.393] 

10.964 

S 

CD 

O 

Velocity  | 

(m/s) 

2.360| 

2.682| 

2.845| 

2.947| 

3.054] 

3.154| 

3.294] 

3.409] 

3.487] 

3.519] 

3.541] 

3.549] 

3.549] 

3.551] 

3.546] 

3.521] 

3.493 

3.415 

3.296 

3.102 

CM 

O 

CO 

Tunnel  3L 

Position  C 

Q398T35H20R1 

c 

-Q 

L_ 

D 

h- 

cC 

17.589| 

13.138| 
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